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Abstract 

Present data on the partonic content of the photon are reviewed. The results 
on the unpolarized structure functions from DIS experiments and on large px 
jet production processes in 77 and 7^ collisions are discussed for both real and 
virtual photons. A few related topics like the QED structure functions of the 
photon and the structure function of the electron are also shortly discussed. 
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1 Introduction 



The concept of the hadronic (partonic) structure of photon is used in describing high 
energy photon - hadron interactions (for a general discussion see e.g. review articles 
I]]). There are two basic types of inclusive processes where the structure of photon is 
tested in existing experiments: 

• the deep inelastic scattering (DIS 7 ), ej — > e hadrons, where the structure func- 
tions of photon are measured, 

• the large pt jet production in jp and 77 collisions (one of the so called resolved 
photon processes), where individual quark and gluon densities in the photon may 
be probed. 

The early experiments of the first type were performed at PETRA and PEP e + e~ 
colliders || [J Final analyses of LEP 1 data on F2 taken at CM energy ~ Mz appeared 
recently, new data are being collected and analyzed at e + e~ LEP collider running at 
higher energies (CM energy: 130-136 and 161-172 GeV), and at KEK collider (energy 
~ 60 GeV). All together they cover a wide range of Q 2 from 0.24 to 390 GeV 2 ; the 
range of the xsj variable extends from ~ 0.002 to 0.98. 

The latter type of measurements, i.e. measurement of the jets in the resolved 
photon processes, started a decade later Q Data have been and are still being taken 
in photoproduction processes at the ep collider HERA (with the CM energy 300 GeV) 
and in 77 collisions at the above mentioned e + e~ machines. They have just started to 
give first results on the light quark (the effective parton density) and gluonic content 
of the photon. 

In this paper the status of recent measurements of the structure function of unpo- 
larized photons for both real and virtual photons in "DIS 7 " experiments and in large 
Pt jet processes involving resolved photons is discussed []. In the presentation we focus 
on recent results (data collected in years 1990 to 1996)0, where qualitative change has 
appeared in both types of measurements. Final results based on few years' runs at 
LEP 1 are being published and a summing - up can be done. On the other hand, 
impressive progress has been obtained in pinning down the individual parton contri- 
butions in resolved photon processes with large pt jets at HERA. These two methods 
of studying the "structure" of photon are becoming recently even more closely related, 
as it turned out that in the reconstruction of K7 by the unfolding procedure the de- 
tailed analysis of the final hadronic state, with the contribution due to resolved photon 
processes, has to be performed. 

The discrepancies observed by many collaborations in the description of the hadronic 
final states in the DIS 7 experiments have enforced the advanced study of various as- 
pects of hadron production in this experimental setup. We found it important to 

1 The first measurement was done in 1981 by the PLUTO Collaboration 

2 In principle the measurement could be performed also at the SLC collider (at energy ~ 91 GeV). 

3 First evidence for the resolved process was found in 1990 by the AMY collaboration at KEK j|. 

4 So far there is no data for the polarized photon; for recent discussion of the polarized structure 
functions of photons see e.g. [|[ ||. 

5 The exception is made for the structure of the virtual photon and the electromagnetic structure 
function for the photon, where old experimental data are also reviewed. 
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collect these results together and in close relation to the large pr jet production re- 
sults, where problems with the proper description of the underlying events appear as 
well. We hope that this way of presentation may help to clarify the situation with the 
production of hadrons in photon induced processes. 

A very comprehensive paper "A compilation of data on two-photon reactions lead- 
ing to hadron final states" by D. Morgan, M. R. Pennington and M. R. Whalley || 
contains data published up to March '94. Beside the inclusive hadron production it 
contains also the data on the exclusive processes and on the total cross sections, which 
are beyond the scope of the present paper. 

Our survey is based on the published data, with few very recent preliminary re- 
sults presented at the 1997 conferences. The short descriptions of results obtained 
by the experimental group together with some representative figures as well as com- 
ments/conclusions quoted from original publications are given. We describe the data 
by the name of the collaboration in alphabetic order, the publication date and the 
reference listed at the end of the paper. Each of the different topics is introduced by 
a short review of the basic theoretical ideas and notation. 

In Sec. 2, data for the real photon are discussed. Sec. 2.1 deals with DIS 7 results 
(structure function F%), Sec. 2. 2 with data on hadronic final states in DIS 7 experiments 
for the real photon and Sec. 2. 3 with large p? jet production in 77 and 7p collisions. 
In Sec. 3 existing data for the virtual photon are presented. In Sec. 4 related topics are 
shortly discussed, namely leptonic structure functions of the photon (Sec. 4.1) and the 
structure function of the electron (Sec. 4. 2). Sec. 5 contains the summary and outlook. 
In the Appendix, the existing parton parametrizations are listed. 
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2 Partonic content of the real photon 



2.1 DIS 7 experiments 

In this section we consider the standard DIS 7 measurements in unpolarized ej collisions 
based on the process (Fig. |1|): 

e7 — > e hadrons, (1) 

at e + e~ colliders (so called single - tagged events) |1J. Here the target is a real (in 
practice almost real) photon, coming from the initial electron (positron). 

e 




X 



Y(pio) 

Figure 1: Deep inelastic scattering e~f (P 2 « 0) — ► eX. The photon probe has virtuality 
q 2 = -Q 2 - 

The flux of these real photons in present experimental setups can be approximated 
by the Weizsacker - Williams formula (called also Equivalent Photon Approximation 
EPA) [|7|, [I]]. Q It corresponds to the typical soft bremsstrahlung spectrum. [] There 
is a possibility to introduce here the structure function of the electron, which would 
contain the convolution of the Weizsacker - Williams flux for virtual photons and the 
structure function of 7* (see discussion in Sec. 3 and 4). 

The cross section for process (|l|), where the real photon with four momentum p 
is probed by the virtual photon with four momentum q = k — k' (q 2 = —Q 2 < 0), is 
given by: 

J -[(l-y)FZ + x Bj y 2 F?\ (2) 



dxBjdy Q A 

'[(i + (i-ym-ft 7 ], (3) 



2ira 2 (2p ■ k) rn M ^, w _,. , 



Q 4 

with the relation among the transverse Fj> (= F?) and longitudinal Fl structure 
functions: F% = F2 + IxBjF^. k(k') are four momenta of initial (final) electron and 

6 Recent critical discussion of this approach for the 77 processes in e + e~ collisions at LEP 2 can 
be founding. 

7 Note that at or 77 options planned at NL Colliders, where the real photons would be obtained 
in the Compton backscattering process, the corresponding flux is expected to be much harder H. 
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the standard DIS variables are defined by 

Q 2 p-q m 

2p ■ q p ■ k 

Note that in the LAB frame y describes the scaled energy of the exchanged photon. 
In practice, the variable y is small and the cross section (2) is effectively saturated by 
F%. Note that in e + e~ collisions there is always a small off-shellness of the photon- 
target, p 2 = —P 2 7^ 0, and there may appear in addition the third structure function, 
F% . It disappears however from the cross section in the single-tag measurements after 
integration over the azimuthal angle (since it enters the cross section as the term 
F3 cos(20)). See Sec. 4.1 for the discussion on the additional structure functions. 

2.1.1 Theoretical description 

In the Quark Parton Model (QPM or PM) one assumes that the hadronic final state 
in Eq. ([I]) is due to the production of quark pairs: q^ and qi (i— 1,2... Nf - number of 
flavours) with the fractional charge Qi. The F£ is obtained by the integration over 
the transverse momentum of the outgoing quark with respect to the target photon 
direction. The full (PM or the lowest order) expression for F2, keeping the terms with 
quark mass m qi , is given by the Bethe-Heitler formula Q: 

a Nf Am 2 
Fl = -N c J2Qtx B ,[(-l + Sx Bj (l - x Bi ) ~ x Bj {l - x Bj )—^)(5 
71 i=i W 



+ [x% + (1 - x Bj f + x Bj (l - 3x Bj )^ - x%^] In i±|], (5) 
where the quark velocity (3 is given by 



Am 2 , 

Hi 



Ami 



XBj 



s \ Q 2 (l-x Bj ] 



(6) 



The energy of the 7*7 collision, s, is equal to the the square of the invariant mass of 
the hadronic system W 2 and is given by 

s = W 2 = ^(l-x Bj ). (7) 

XBj 

In the limit of s well above threshold, i.e. for /3 ~ 1, which for fixed x B j (not too 
small and not too close to 1) corresponds to the Bjorken limit, one gets 

1-13 1-P mlxBj 
The structure function F2 can be approximated in such case by 

a Nf W 2 
F 2 = -N e E Qi x Bj[[x Bj + (1 - x Bj ) 2 } In — + 8x Bj {l - x B j) ~ 1] (9) 
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and it can be used to define the quark densities in the photon: 



2N f 

n = XB J EQhUxB 3 ,Q 2 ). (10) 
i=l 

In the above formulae N c =3 denotes the number of colors. In the last equation a 
(natural) assumption, that quark and antiquark distributions in photon are the same, 
has been introduced. Note that 

• F2 is calculable in the PM, in contrast to the structure function of hadrons, e.g. 
the nucleon structure function F 2 N . 

• F2 is proportional to a, the fine structure coupling constant (« 1/137). There 
is an overall logarithmic dependence on the energy scale squared W 2 (or Q 2 , 
see Eq. (^)). A large value of F% , or in other words a large quark density, is 
predicted at large x B j. 

• F2 is a sum of the quark densities with the factors Qf , where Qi is fractional 
charge of quark, so the individual quark density is proportional to Q 2 (see Eqs. 

(B W- 

The longitudinal structure function F2 is not zero in PM, in contrast to the corre- 
sponding function for hadrons, and it is scale invariant: 

N f 

Fl = -N e Y,Qtx Bj [4x Bj (l -x Bj )}. (11) 
71 i=i 

Note that also the third structure function is scale invariant: 

N f 

ry J 

= -N c ^Qix Bj [-x 2 Bj }. (12) 
n i=i 

In the leading logarithmic approximation (LLA) the Callan-Gross relation Fq, = 
2x B jFi holds as in the case of hadrons. 

In this approximation the PM formula for the quark density is given by 

Q7(xBj, Q 2 )\ L p L m = ^N c Q 2 [[x 2 Bj + (1 - x Bj ) 2 } In -^H. (13) 

Since the LL contribution corresponds to the on-shell quarks one can treat x B j in the 
above formula as equal to the part of four momentum of the initial photon carried 
by the quark; this latter variable is usually denoted by x 1 (or simply x). Note also 
that in the above formula instead of a quark mass there appears the QCD scale Aqcd- 
Therefore this expression allows to describe all the light quark (light as compared 
to the scale W 2 ) contributions to F2 in a universal way. Heavy quark contributions 
should be treated separately, according to the Eqs. (^|, ||) [[|. 

An additional Q 2 dependence will appear in and in q] due to the QCD correc- 
tions, which can be described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equa- 
tions or by other techniques [I], |TD|, [TTJ, [T2], |T3| . In this framework the gluonic content 
of photon, g 7 , appears as well. 
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The inhomogeneous DGLAP equations for the real photon can be represented in 
the following way (below q] is used for quarks and antiquarks): 

W" + §r d J lP " {X y ) I qUy) + P " {X y )g " (y)] (14) 

with the standard splitting functions P qq , P qg , P gq , P gg and in addition with the function 

P qi = N c [x 2 + (l-x)% (15) 

describing the splitting of the photon into quarks. 

Note that it is possible to solve the above equations without the initial conditions 
TO] , assuming the LL or NLL behaviour of the solution. Obtained in this way the, 



so called, asymptotic solutions, have a singular behaviour at small x [13[]. Therefore 
in practice while solving the equations (14) the initial conditions have to be assumed 
from a model or taken from measurements at some (low) Ql scale ||14j| . On this basis 



the parton parametrizations for the photon can be constructed, see the Appendix for 
the details of existing parton parametrizations. 

In pre-QCD times, the hadronic structure of the photon has been solely attributed 
to the vector meson component (p, in the (real) photon. Therefore for the matrix 
element of the photon between e.g. the nucleon states, the following representation 
by the corresponding matrix elements for the p current was assumed ||15|: 



m 2 1 

<N\J IM \N>= p — < N | Plt \ N > +... (16) 

9 P q - m* 

That was the basic assumption of the Vector Dominance Model (VDM) or Gen- 
eralized VDM (GVDM) if higher vector mesons states were included [|T^, [IJ|. This 
non-perturbative (or hadronic) contribution to the structure of the photon is present 
of course in the measured structure function as well as in q 1 or g 1 extracted from 
the experimental data. 

2.1.2 Data on F 2 7 

Here we discuss the latest results for for a real photon based on single-tagged events 
at e + e~ colliders {i.e. DIS 7 events, see Eq. (U) and Fig. p. Note that in practice not 
XBj but the quantity x vis is measured, 

Xvis = Q 2 + wi; (17) 

where W V i S is the invariant mass of the visible hadronic system. For the nonzero 
target-photon virtuality P 2 , the measured quantity is 

Q 2 



X i ; 



Q2 + W 2 s+ p2- 
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From a measured value for x V i S one can reconstruct xsj by an unfolding procedure. 

Recently, in order to perform reliable unfolding, the events are grouped into classes 
of similar topology (see e.g. DELPHI 96b), which roughly coincide with the con- 
tributions due to QPM, VDM or RPC (RPC - Resolved Photon Contribution, where 
photons interact through their partonic contents). 

In the early analyses the separation between these contributions (QPM with QCD 
corrections and VDM, basically) was obtained by introducing the parameter of the 
minimal transverse momentum of the produced hadrons, Pj>- The VDM (or hadronic) 
contribution should be negligible for the production of particles with p? larger than pj,, 
where, on the other hand, other contributions should be included. The FKP approach 



17 1, where this kind of cutoff parameter is built in, was natural here. Note that 



in generation of two quark final state (in e.g. F2GEN event generators) two extreme 
types of angular dependence for final quarks are assumed: point-like (PL) with angular 
dependence as in the lepton pair production in real photon-photon collision and a 
"peripheral" one with the distribution corresponding to an exponential form of the 
Pt dependence, with a mean of 300 MeV (as if the photon interacted as a hadron) []. 
Nowadays the Monte Carlo generators HERWIG and PYTHIA, adapted for DIS 7 in 
1995, can describe all these types of contributions with any of the existing parton 
parametrizations. It is no longer necessary to fit an empirical p 1 ^ parameter to the 
data before unfolding. 

It is worth noticing here, that at present energies (or Q 2 scales) light quark (it, d, s) 
and heavy quark distributions are treated differently. For c, b quarks the QPM formula 
(Eq. (^j)) is applied (at least close to the threshold), while for the lighter quarks the 
QCD corrections are necessary. Note also that in the earlier measurements the c quark 
contribution was usually subtracted from the structure function F£. 

Unfolding can be performed using the traditional method based on linear scale or 
using the new approach with the logarithmic scale, especially useful for extracting F7 
results in the small xsj region. 

As we have already pointed out, the study of the hadronic final state became part 
of the measurements of the structure function for the photon at e + e~ colliders. The 
details of the studies of the hadronic final state will be presented in the next section. 
Here we would only like to stress that some discrepancies have been found for certain 
distributions, like the pseudorapidity r\ of final hadrons distributions. This fact is 
included in the estimation of the uncertainty of the measured function F^. 

The general features of FJ, as far as Xsj and logQ 2 dependences are concerned, 
agree with the theoretical expectations, although the precision of the data does not 
allow in many cases to distinguish between existing parton parametrizations and/or 
clarify the small XBj behaviour of F%. 

We start the presentation of the data from the LEP collider, then TRISTAN data 
are discussed. Collective figures of versus xgj and F7 versus Q 2 are presented at 
the end of this section. 

The sets of data are described by the name of the collaboration (in alphabetic or- 
der), the publication date, and reference listed at the end of the paper. In comments 
we quote statements from the original papers (for abbreviations used for the parton 

8 also a mixed 'perimiss' dependence was studied |p3. 
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parametrizations see the Appendix). 



DATA 



•ALEPH 97a [L8j (LEP 1) 

Data on r 7 ^ at LEP 1 energy for xsj from 0.002 to 0.9 and Q 2 between 6 and 44 GeV 2 
were collected in the period 1991-94. An analysis of the hadronic final state was per- 
formed using the QPM+VMD model and the HERWIG simulations (with GRV, GS 
and LAC1 parton parametrizations) (see next section). 



ALEPH <(f>=!L9 GevVc 


PRELIMINARY 




, GRV 




LAC1 




MID 





ALEPH <Q->=I9.1 GeTIc 


- | PRELIMINARY 








GRV 




LAC1 




SaSlD 





0,5 0.6 0,7 



Figure 2: The unfolded F% /a for < Q 2 >=8.9 and 19.1 GeV 2 compared with predictions 
of GRV, LAC1, and SaSlD parametrizations. 

The unfolded results for F2 are presented in Fig. 0. F2 averaged over xsj and Q 2 
in two bins of Q 2 is given in the table: 



<Q 2 > 
[GeV 2 ] 


x Bj 


(stat. + syst.) 


8.9 


0.3-0.8 


0.36 ±0.16 ±0.06 


19.1 


0.3-0.8 


0.44 ±0.08 ±0.02 



Comment: Problem with the proper description of hadronic energy flow as a function 
of pseudorapidity or azimuthal separation angle was found (see next section). 

•ALEPH 97b |]| (LEP I) 

Data taken from 1991 to 1995 at average Q 2 equal to 279 GeV 2 . The hadronic final 
state was also studied using four QPM±VDM models and HERWIG Monte Carlo 
program (with GS2, GRV LO, LAC1 and SaSld parton parametrizations, see next 
section). The results for F2 are presented in Figs. |3] and 

Comment: Problem with the proper description of the pseudorapidity distribution of 
final hadrons (see next section). 
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x 



Figure 3: The structure function F^ /a as a function of XBj (from 




Q 2 (GeV 2 /c 2 ) 



Figure 4: The Q dependence of F^/a for the range 0.3< XBj <0.8 (from JT^/ j 



•DELPHI 96a @ (LEP 1) 

Data on F2 were taken in the period 1991-93 for Q 2 between 4 and 30 GeV 2 and 
for %Bj down to 0.003. The so called F] D ^ was also measured and compared 
with QED prediction, also for nonzero averaged virtuality P 2 of the target photon 
(see Sec. 4.1 for details). Estimated target photon virtuality P 2 was used in the 
unfolding of i 7 ^ Q The TWOGAM event generator was used to simulate QPM events 
and another event generator was used to obtained the QCD correction (LL) to the 
point-like contribution for light quarks in the FKP approach. The GVDM and the 
point-like (FKP) contributions were studied, with =0.1 and 0.5 GeV. Results for 

9 It was found that although < P 2 >=0.13 GeV 2 , a fixed value of 0.04 GeV 2 fits the data better. 
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Fl at < Q 2 >=12 GeV 2 are presented in Fig. | and in the table below: 



<Q 2 > 
[GeV 2 ] 


X Bj 


Fq/a 
(stat. + syst.) 


12 


0.003 - 0.080 
0.080 - 0.213 
0.213 - 0.428 
0.428 - 0.847 


0.21 ±0.03 ±0.06 
0.41 ±0.04 ±0.05 
0.45 ±0.05 ±0.05 
0.45 ±0.11 ±0.10 





0.8 






a 




a 






0.7 


a 














0.6 








0.5 




0.4 




0.3 




0.2 




0.1 








DELPHI <Q 2 > = 12 GeV V 



OPAL <Q> = 14.7 GeV 7c 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.i 



Figure 5: Unfolded F^/a for the light quarks (DELPHI and OPAL data). The curves show 
the sum of the GVMD model prediction multiplied by the threshold factor 1-XBj and the 
prediction of the FKP parametrization for the point-like part of F^, with different values of 
parameter p^: 0.1 GeV (upper line) and 0.5 GeV (lower line). The bottom curves show the 
GVDM contribution with different target masses (from |2(]/). 



The averaged value of F^jcn over the %Bj range between 0.3 and 0.8 was extracted: 



<Q 2 > 
[GeV 2 ] 


< Fq/a > 


12 


0.45 ±0.08 



For comparison with other measurements of the Q 2 dependence of the F% , see Fig.|6]. 
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2.5 



a 

S 2.25 



^ 1.75 



1.5 
1.25 
1 

0.75 
0.5 
0.25 




□ PLUTO 

* JADE 

A TPCKl 
t TASSO 
a AMY 
> OPAL 

• DELPHI 



0* 



0.3<x<0.8 

tP t = 0.1 GeV/c 

= 0.27 GeV/c 
\l\ =0.5 GeV/c 
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10 , , ^° 

g 2 (GevV/j 



Figure 6: F^/a averaged over XBj between 0.3 and 0.8. The curves show the FKP 
parametrization predictions for different values of the parameter p\ (from f&(\]). 



Study of the behaviour at < Q 2 >=12 GeV 2 in the low XBj domain leads to 
following results (Fig. |7|): 



<Q 2 > 
[GeV 2 ] 


x Bj 


Fq/* 

(stat. + syst.) 


12 


0.003 - 0.046 
0.046 - 0.117 
0.117-0.350 


0.24 ±0.03 ±0.07 
0.41 ±0.05 ±0.08 
0.46 ±0.17 ±0.09 



DELPHI 

-1 i 




i 




GRV 


\ 




GS 






DO 




LAC1 





0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 7: Unfolded F% /a for < Q 2 >= 12 GeV 2 from the DELPHI experiment with the 
LO parametrizations LAC1, GS, DO, and GRV (from fl2fj/). 



Comment: No rise of F% at small %Bj has been found. GRV and GS leading order 
parametrizations of the quark density in the photon are in agreement with data. 



•DELPHI 96b J21| (LEP 1) 



Measurement of the photon structure function (data collected in the years 1991- 
95) together with a study of the hadronic final state, in particular the resolved target 
photon contribution (see next section), is reported. 
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Two types of models for the extraction of the VDM contribution were used for the 
< Q 2 >=13 GeV 2 data (GVDM and "TPC/2 7 "-type, for details see Ref.0). QPM, 



and RPC contribution with the GS2, SaS4 and GRV3 parton parametrizations were 
studied. The results for the < Q 2 >=106 GeV 2 were also obtained. The unfolding; 
was done in a linear and logarithmic scale in XBj- 

Fig. |8| shows the xbj dependence for both Q 2 samples. 





Figure 8: Unfolded F^/a compared with QPM+GVDM+RPC (GS2) (solid line) and pre- 
dictions of SaS4 (dashed line), GS2 (dashed-dotted) and GRV3 (dotted) parametrizations 
(from JH/J. 

The averaged value of Fj /a. for 0.3 < xbj < 0.8 as a function of Q 2 was extracted: 



<Q 2 > 
[GeV 2 ] 


< Fj/a > 
(stat. + syst.) 


13 
106 


0.38 ±0.031 ±0.016 
0.576 ±0.081 ±0.076 



Comment: The importance of the final hadronic state topology was noticed and the study 
of a linear and logarithmic unfolding performed. 

•DELPHI 97 fl2|] (LEP 1, LEP 2) 

This paper reports a recent study of in the Q 2 range between 3 and 150 GeV 2 , 
based on data from the 1994-95 runs for energies around the Z° mass and from 1996 
for energies between 161 and 172 GeV. An analysis of the hadronic final state is per- 
formed and compared with predictions of the TWOGAM Monte Carlo program, where 
QPM, VDM and RPC(GS2) parts are included (see next section). 

Unfolded results for F^ as a function of xsj and Q 2 are presented in Figs. [|a 
and [|b, respectively. 

Comment: The data for hadronic final state distributions, including energy flow versus 
pseudorapidity, agree with predictions of the TWOGAM simulations. 
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0.4 
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LEP1 and LEP2 preliminary b) 

o LEP1 
T LEP2 



LEP1 

Previous analyst 



1.2 1.4 1.6 1.8 



lo g 10(<?) 

Figure 9: a) F^/oc versus XBj from the DELPHI experiment based on LEP 1 data 
for three values of < Q 2 >: 6.5, 13, and 22 GeV 2 . The solid line corresponds to 
QPM+GVDM+RPC(GS2), dotted - QPM+GVDM and dashed - GVDM. b) F^/a aver- 
aged over XBj from 0.3 to 0.8. Results obtained from LEP1 and LEP2 data are shown 
together with results of a previous analysis of LEP1 data (from /p3/). 



•OPAL 94 [H (LEP 1) 

Measurement of at < Q 2 >=5.9 and 14.7 GeV 2 was performed using data from 
the period 1990-92. The VDM and pointlike contributions (in the FKP approach) are 
separated by the cutoff parameter pj, found to be 0.27±0.10 GeV. 



Results compared with PLUTO and TPC/27 data are presented in Fig. |10 



OPAL <Q"> = 5.9 GeV 



□ PLUTO <Q 2 > = 4.3 GeV ! 



A TPC/27 


4 < Q* < 7 


GeV' 










- I 
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• OPAL <0 2 > = 14.7 GeV ! 
□ PLUTO <0"> = 9.2 GeV 2 



0.3 G.4 0.5 0.5 0.7 




0.1 0.2 0.3 0.4 0.5 0.5 0.7 0.8 0.9 



Figure 10: Unfolded F^/a shown with previous measurements in other experiments at 
similar < Q 2 > values. The curves show contributions of VDM (dots), QCD-based model 
(dashes) and their sum (line) for a) < Q 2 >=5.9 GeV 2 b) < Q 2 >=14.7 GeV 2 (from pfj. 



•OPAL 97a HH (LEP 1) 

The measurement of F 2 7 was done for 6 < Q 2 < 30 GeV 2 and 60 < Q 2 < 400 GeV 2 , 
using the full sample of data at the Z peak (years 1990-95). The detailed analysis of the 
hadronic final states was performed and sizeable discrepancies with the expectations 
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were found especially at low xsj (see next section). The influence of the choice of 
different MC generators (HERW1G, PYTH1A) on K7 was studied. The estimation of 
the P 2 values was done based on the SaSl parametrization. 

The results for the F^, presented in Fig. |TT|a,b,c, are as follows (the value of F% /at 
is given at the centre of the XBj bin): 



<Q 2 > 
[GeV 2 ] 


x Bj 


Fq/a 
(stat. + syst.) 


7.5 


0.001 -0.091 
0.091 -0.283 
0.283 - 0.649 


0.28 ± 0.021&5S 
0.32 ± 0.02±8;?1 
0.38 ± 


14.7 


0.006 -0.137 
0.137-0.324 
0.324 -0.522 
0.522 -0.836 


0.38 ±0.011^3 
0.41 ± 0.02l8:oi 
0.41 ± 0.03±g;SS 
0.54 ±0.05lg:?3 


135.0 


0.100 -0.300 
0.300 -0.600 
0.600 -0.800 


0.65 ± 0.091^ 
0.73 ± 0.08±8;g| 

0.72 ± eloign 
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Figure 11: F^/ol data from the OPAL experiment for the number of flavours Nf=4. Curves 
in (a)-(d) show predictions of GRV and SaS ID parametrizations (from 



The measurement of F^/a as a function of Q 2 averaged over x B j range, 0.1< 
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XBj <0.6, leads to following results (Fig. |TT|d) : 



<Q 2 > 
[GeV 2 ] 


< Fq/a > 
(stat. + syst.) 


7.5 
14.7 
135.0 


0.36 ± 0.02i^ 
0.41 ± 0.01±g;8| 
0.71 ± 0.07ig;Jf 



and the slope d{F^ /a)/dlnQ 2 is measured to be 0.13lo;o4. 

Comment: No correction for P 2 ^ was made. Large discrepancies between the 
hadronic energy flow data and Monte Carlo simulations are observed at low x V i S! when 
the results are presented versus pseudorapidity or azimuthal angle (see next section). 



•OPAL 97b (LEP 1) 

Measurement of F 2 7 at LEP 1 was done for < Q 2 >=1.86 GeV 2 and < Q 2 >=3.76 
GeV 2 as a function of XBj, reaching the lowest ever measured (center of bin) value: 
££j=0.0025. For a better sensitivity on the low XBj region the unfolding procedure was 
performed on a logarithmic scale. Final state topology was analysed as well (using the 
HERWIG, PYTHIA, F2GEN generators - both with the pointlike and the peripheral 
distributions). Discrepancy between the data (hadron energy flow) and results from 
the Monte Carlo generators, as well as between different models was found for low 
x vis <0.05 (see Sec. 2. 2). 

Obtained values of F% ja are shown in Fig. [12] together with measurements from 
PLUTO and TPC/2 7 . 



Comment: "No correction for P 2 ^0 was made. GRV-HO is consistent with the low-x 
OPAL results in the lower Q 2 bin, but at higher Q 2 it underestimates the low x OPAL 
data. GRV-LO and SaSID describe the unfolded results worse than GRV-HO. Shapes 
of measured are flat within the errors, but a small rise in the low x region is not 
excluded. " 
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Figure 12: OPAL F% /a data (circles) as a function of XBj for < Q 2 >= 1.86 GeV 2 (a) and 
< Q 2 >= 3.76 GeV 2 (b). Also data from PLUTO (crosses) for < Q 2 >= 2.4, 4.3 GeV 2 and 
TPC/2 r y (squares) for < Q 2 >=1.31, 2.83 GeV 2 are shown. The curves show predictions of 
GRV-HO (dots), GRV-LO (line), and SaSID (dashed) parametrizations. The range of the 
XBj-bins of the OPAL results are marked at the tops of hgures (from figjjj). 



: i i i 

• OPAL, <Q 2 ) = 3.76 GeV 2 



o PLUTO, ( Q 2 > = 4.3 GeV 2 




•OPAL 97c pH (LEP 2) 

New data on from LEP 2 at the CM energies 161-172 GeV were collected in 1996. 
The < Q 2 > range lies between 9 and 59 GeV 2 . Also distribution of the final hadronic 
energy flow was studied (HERWIG, PYTHIA and F2GEN, see next section). 

The unfolded /ot as a function of XBj and Q 2 are presented in Figs [13], [L4]. A 
special study of the Q 2 dependence of is performed for different XBj ranges, see 



Fig. 15 



A fit to the new data at the energies 161-172 GeV and the previous OPAL set 
at 91 GeV for Q 2 from 7.5 to 135 GeV 2 , averaged over the XBj range of 0.1-0.6 (see 
Fig. [14]), has the form: 



F 2 7 (Q 2 )/« = (0.16 ± 0.05lg;$ + (0.10 ± 0.02^) ln(Q 2 /GeV 2 



Comment: No correction for P 2 ^0 was made. Discrepancies are observed in the 
hadronic energy flow between the data and the HERWIG and PYTHIA simulations, 
especially at x vis <0.1. Accuracy of the data does not allow to see the expected different 
slope of F2 versus logQ 2 for different xgj ranges. 



16 



» 1 
So.75 
0.5 
0.25 





GRV 

I i i i I i i i I i i i r v 



J 1 

: <Q 2 > = 30.GeV 2 


(b) 


r OPAL 








< 








ASYM 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

x x 



8 1 



< Q 2 > = 14.5 GeV 2 (c) 



< Q 2 > = 59. GeV 2 



(d) 



-0.75 
0.5 
0.25 




OPAL 



OPAL 



s 1 
®0.75 
0.5 
0.25 


0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

x x 



I ... I ... I i'i'VI- 



i . . . i . . . rr-rv-i- 




J 1 

| <Q 2 > = 41.GeV 2 

r OPAL 


(f) 







0.2 0.4 0.6 0.8 1 

x 



Figure 13: XBj-dependence of F^/ot measured by OPAL for different value of < Q 2 >: 9, 
30, 14.5, 59, 11, and 41 GeV 2 . For comparison the GRV LO (solid line) and SaSID (dashed 
dotted line) parametrizations are shown; dotted line represents the hadronic component 



(VDM?) (HAD) and the dashed one the asymptotic solution (ASYM) (from [26]). 
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Figure 14: F% /a averaged over 0.1< XBj <0.6, for the energy 91 GeV (open circles) and 
for the energies 161, 172 GeV (full circles). Predictions of the QCD calculation are shown 
by the lines: solid (GRV LO), dotted (SaSID), and double-dotted (HO, based on GRV 
HO parametrization for light quarks); dashed line corresponds to the asymptotic solution 
(ASYM). The charm contribution is calculated separately (from pq/J. 
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•AMY 95 |27f (TRISTAN) 

A high Q 2 measurement of the photon structure function K7 was performed. Unfolded 
F2 is presented in Fig. [16] (below are the corresponding numbers). 



<Q 2 > 


< X B j > 


FZ/a 


[GeV 2 ] 




(stat. + syst.) 


73 


0.25 


0.65 ±0.08 ±0.06 




0.50 


0.60 ±0.16 ±0.03 




0.75 


0.65 ±0.11 ±0.08 


390 


0.31 


0.94 ±0.23 ±0.10 




0.69 


0.82 ±0.16 ±0.11 



The values of averaged over 0.3< xbj <0.8 are as follows: 



<Q 2 > 
[GeV 2 ] 


< nioi > 


73 
135 


0.63 ±0.07 
0.85 ±0.18 



Comment: The observed xbj -behaviour of F2 is consistent with the GRV parametriza- 
tion and with the FKP one (fitted parameter pj, is equal to 0.51±0.39 GeV, if only 
AMY data are included, and 0.45±0.07 GeV if all available, at that time, F^ data are 
taken). 
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Figure 16: F^ /a data from the AMY collaboration compared with a few parametrizations 
of the photon structure: FKP, QPM (quark parton model), VMD model fab), WHIT1, 
LAC1, DG fad), GRV, AFG faf). Upper fac,e) and lower (b,d,f) figures correspond to 
averaged Q 2 =73 GeV 2 and Q 2 =390 GeV 2 , respectively (from HJ/J. 



•AMY 97 (TRISTAN) 

The measurement of F% at < Q 2 >=6.8 GeV 2 was performed; the results for F^ versus 
Xsji together with earlier data, are presented in Fig. [L7|. The comparison with various 
parametrizations is given in Fig. |18| . 



K PLUTO Q ! - 4.3 GeV ! □ TOPAZ Q ! = 5.1 OeV* 

o PLUTO Q* - 8.2 GeV 1 X OPAL q" = 6.9 GeV* 

•> TPC/27 Q ! - S.l GeV 1 • AMY (j' - 6.8 GeV" 




Figure 17: The measured F% /a values compared with results of other experiments at Q 2 
4-9 GeV 2 (from |2|/j. 
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Figure 19: TOPAZ 94 data for /a versus XBj for Q 2 =5.1, 16, 80 GeV 2 and a comparison 
with parton parametrizations: a) FKP, QPM (parton model), VMD model and b) DG, 
LAC1, GRV (from $MI). 



The results are presented in Fig. [19]. To study the Q 2 dependence of K/, the 
averaged values of F^/ot in the a^-range from 0.3 to 0.8 were extracted at < Q 2 >= 
16, 80 and 338 GeV 2 (the numbers in parentheses are the results for light quarks 
alone): 



<Q 2 > 


Fq/a 


[GeV 2 ] 




16 


0.47 ±0.08 (0.38 ±0.08) 


80 


0.70 ±0.15 (0.49 ±0.15) 


338 


1.07 ±0.37 (0.72 ±0.37) 



Comment: The final hadronic state described by QPM(c) + VDM+FKP(u,d,s) with 
p%=0.l,0.5 and 1 GeV, also by GRV, DG and LAC1 parton parametrizations, was 
studied. 



For an overall comparison the collective figures of F^jct versus XBj (Fig. ^ p0[ ) 
and F^/at versus Q 2 (Figs. |2l] fl3i"l , |22| [|32]1 ), containing also earlier data not discussed 
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here, are presented. In Fig. |20| the comparison with theoretical predictions for in- 
dependence of F2 based on SaS-lD (LO) and GRV NLO parton parametrizations are 
shown (note the logarithmic scale). For comparison with theoretical predictions of the 
Q 2 evolution for different xsj ranges see Fig. |22], and also Figs. f|, |6], |9|b, |TT|d, |T3| , [L5. 

The effective parton density, as measured at HERA collider by HI group (see HI 
97a |33[) in the jet production from resolved photons, is compared with the data 



in Fig. 23 as well. 
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Figure 20: The photon structure function ¥^ /a as a function of x Bj in bins of Q 2 compared 
to the GRV NLO (solid line) and SaS-lD (LO) (dashed line) parametrizations of parton 
distributions in the photon (from pfl/J. 
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Figure 21: Tie photon structure function t^/ot (data from various experiments) as a 
function of Q 2 (from Mj). 
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Figure 22: The Q 2 evolution ofF^/a measured by OPAL (0.1 < x B j < 0.6), AMY (0.3 < 
x B j < 0.8), JADE (0.1 < x B j < !■), DELPHI (02 < x B j < 0.8), TPC (0.3 < x B j < 0.6), 
TOPAZ (0.3 < xsj < 0.8). The lines are predictions of parametrizations averaged over 
different XBj ranges, HO calculation being a sum of GRV HO parametrization for light 
quarks and the charm contribution calculated independently from GRV; also the asymptotic 
solution is shown (ASYM, dashed line) (from ft3$H). 
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Figure 23: The scaling violation of the parton densities compared with the Q 2 dependence 
of the F2 (averaged over the range 0.3< xsj <0.8) (from 
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2.2 Hadronic final states in the DIS 7 experiments 

Although a detailed analysis of hadronic final states, including those involving large pt 
particles/jets, has not been the main aim of DIS 7 experiments, nowadays it has proved 
crucial in extracting the unfolded F%. The obtained experimental results, and problems 
that appear while describing hadronic energy flow, transverse energy, pseudorapidity 
and other distributions within QCD, deserve close attention and a separate treatment 
- this section is devoted to this subject. 

It is worth noticing that as corresponds to the total cross section for hadron 
production in 7*7 collisions, all the regions: Q 2 <C p?, Q 2 ~ Pt an d Q 2 ^ Pt 
contribute, with the bulk of the cross section being due to the soft production, i.e. 
with a relatively small pt- In the next section dedicated analyses of the large pt jet 
production in real 77 and jp processes will be presented. Some of the problems that 
will appear there are common with these presented here. 



2.2.1 Modeling of the hadronic final state in 7*7 collision 

In this section the 7*7 collisions resulting in the hadronic state will be discussed. As 
we have already mentioned in Sec. 2. 1.2, the QPM, VDM and RPC contributions are 
introduced to model the hadronic final state. In recent analyses the MC generators 
HERWIG and PYTHIA are being used with chosen parton parametrizations. In the 
OPAL analysis, the generator F2GEN (based on TWOGEN used in older analyses), 
with assumed two quarks in the final state is used in addition, for the comparison and 
systematic checks. In this generator one can easily put a specific assumption on the 
final states. For example different angular distributions may be assumed in the 7*7 CM 
system: pointlike (PL), i. e. as for lepton pair from two real photons, peripheral with an 
exponential distribution of the transverse momentum, or the "perimiss" combination 
II. In a recent DELPHI analysis, the TWOGAM generator ("fixed" in 1993) was 



successfully used, with QPM, VDM and RPC contributions included [^1], [2! 

As far as data on global hadronic variables like W V i s , Q 2 , Mij (invariant mass of 
two jets) and other distributions are concerned, there is a fairly good description of the 
data by existing MC generators. The problems arise for the transverse energy out of 
the tag plane (the plane defined by the initial and tagged electrons) or for the energy 
flow per event as a function of pseudorapidity, rj = — ln(tan(6*/2)), where the polar 
angle is measured from the direction of the target photon. The discrepancies are very 
pronounced especially for small XBj- 

The first observation of the disagreement between the data and MC models was 



made by the OPAL collaboration [f2J]. "The serious discrepancies between the data 
and any of the available Monte Carlo models are seen both within the central region 
of the detector (|?7| < 2.3), where the energy flow is well measured, and in the forward 
region, where the energy can only be sampled." It is clear that the unfolding of F% 
"will have large errors as long as the energy flow from different models remains in clear 
disagreement with the energy flow in the data, in particular in the region of x V i S <0.1 
and Q 2 <30 GeV 2 " (from [|4|])- A similar effect has been seen by now by other groups. 



The method used recently to improve the agreement with the data is the so called 
"HERWIG + power law p" generator, where the additional power law spectrum of 
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the form dk\j(k\ + k$) was introduced (it seems to be needed at HERA (ZEUS) as 
well, see next section). 

One should be aware that in hadronic final state in DIS 7 experiments two types of 
large scale may appear: Q 2 and p 2 -,, describing the transverse momentum of the final 
hadrons/jets. The bulk of the data corresponds to events with not very large pt', if 
for these events the relation Q 2 3> p\ 3> P 2 holds, the interpretation in terms of the 
photon interaction between one direct (7*) and one resolved real (7) photon may be 
introduced. Then it is not clear what scale should be used in the parton density for a 
real photon f 1 (x,Q 2 ): pt or Q 2 . Moreover the processes corresponding to Q 2 <C p\ 
should be treated as being resolved from the point of view of both: a real and a virtual 
photon. 

Note, that in recent analyses the dependence on the number of jets in the final 
state is studied. 



2.2.2 Data on the hadronic final state 

The analyses of the hadronic final state accompanying the measurements of dis- 
cussed in Sec. 2. 1.2 as well as the results from independent analyses of large pt hadron 
production in single tagged events are presented below. 

For the two photons involved, we introduce the following notation for their squared 
(positive) virtualities: P 2 = Q 2 and P 2 = P 2 (with P 2 > P|), where the correspond- 
ing variables used in the discussed DIS scattering are also given. 

DATA 



•ALEPH 97a [L8j (LEP 1) 

Hadronic energy flow in the azimuthal angle separation <p sep and pseudorapidity 77 was 
measured for Q 2 between 6 and 44 GeV 2 . The results are presented in Fig. 




Figure 24: Energy flow as a function of 4> sep (a) and rj (b). Predictions: QPM+VDM (solid 
line) and HERWIG generator with GRV (dashed line) (from jdW). 



Comment: "Discrepancy with MC models (QPM + VDM and HERWIG (with parton 
parametrizations GS, GRV, LAC1 ) ) when energy How is plotted as a function of ra- 
pidity and azimuthal separation." 

•ALEPH 97b JT§ (LEP 1) 

Data taken from 1991 to 1995 for the average Q 2 =279 GeV 2 . Four different QPM+VDM 
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models were used and in addition the HERWIG 5.9 generator with the GS, GRV, LAC1 
and SaSld parton parametrizations. 

The pseudorapidity distribution compared with the model predictions is presented 
in Fig. |25|. A disagreement is found, similar to the one observed previously by ALEPH 
at lower Q 2 (ALEPH 97a), but now at the negative pseudorapidity. 




Pseudorapidity 



Figure 25: Pseudorapidity of all charged tracks and neutral calorimeter objects. Histograms 
are predictions of three of the models (dotted line from HERWIG, others - QPM+VDM type) 
(from (TJ/J. 



Comment: "The pseudorapidity distribution can not be described by the models as 
observed previously at lower Q 2 , note however that the excess is now in the opposite 
hemisphere than at lower Q 2 ." 

•ALEPH 97c g (LEP 1) 

Single tagged events collected in the years 1992-1994 at the average Q 2 =14.2 GeV 2 . 
A dedicated study of the final hadronic state assuming the QPM+VDM for u,d and 
c quarks and using HERWIG 5.9 generator (with GRV parametrizations). Also the 
modified HERWIG generator was used which corresponds to the "HERWIG + power 
law p t " , with kto = 0.66 GeV (see also OPAL 97e where this model was introduced). 
Results are presented in Figs. |2B], [27|, The approach "HERWIG + power law p t " 
leads everywhere to better description of the data. 




O 2 4 6 8 10 12 14 16 18 20 
E t Out of the tag plane (GeV) 



Figure 26: The energy (out of plane of the tag and the beam) distribution (from j3^/ j. 

In order to pin down the source of observed discrepancy the analysis of the number 
of cone jets was performed in the final hadronic state (see OPAL 97d for the first 
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analysis of this type). HERWIG model where resolved photon processes are included 
should in principle give larger number of two jet events than QPM+VDM model, 
having additional production mechanisms. On the contrary, it gives less (see Fig. p8|) . 

Also the energy not assigned to jets in two jet events, ^ on ~J e * ; was studied. To 
check the presence of the single and double resolved photon events, x 1 and x tag distri- 
butions were measured. The "HERWIG + power law p t " approach was used sucessfully 
to describe these data as well. 
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Figure 27: The energy Row versus pseudo-rapidity of the final hadrons (from p4H ). 
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Figure 28: (from ffij). 

Comment: " "HERWIG + power law pt" is better in modeling the region of large E t>out 
and the peak at positive rapidity". 
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•DELPHI 95 |3§] (LEP 1) 

First evidence of hard scattering process in the single-tagged ej collision in the data 
from the 1991-1992 run is reported. The values of Et of observed jets were larger than 
1.5 GeV, while the magnitude of the mass of virtual photon was equal to < P 2 >~ 
0.06 GeV 2 . This corresponds to the standard resolved (almost real) photon process, 
from the point of view of both photons. 

Comment: This is not typical DIS experiment since the photon probe is almost real. 
The analysis is not typical, either, for the standard large px jet study. "The data 
are consistent with the predictions for quark and gluon density functions in the GS 
parametrization. The sum of the contributions: VDM + QPM (Quark Parton Model) 
+ [QCD - RPC(Resolved Photon Contribution)] is needed in order to describe the data; 
the DO and LACS parametrizations do not adequately describe the data." 

•DELPHI 96b |2Tfl (LEP 1) 

The results for the averaged Q 2 =13 and 106 GeV 2 were considered and compared 
with TWO GAM generator. The QPM, GVDM contributions and, for the first time 
at DELPHI, RPC with GS2 parton parametrization were included. 

Comment: The resolved photon contributions are needed to obtain the description 
of the data, "...the correct unfolding procedure leads to the pointlike plus hadronic 
part of F% in a shape that is inconsistent with QPM (FKP) + GVDM in the region 

of loW-XBj" 

•DELPHI 97 Q2|] (LEP1, LEP2) 

Analysis of hadronic final state in the K7 measurement (see Sec. 2. 1.2) was performed 
for energies around the Z° mass and for 161-172 GeV (1996 run). The Q 2 ranges: 6, 
13, 22 and 17, 34, 63 GeV 2 were studied, respectively. Hadronic final state topology 
with events containg jets was studied using TWOGAM generator. For results see 



Fig. 29 for LEP2 data and Fig. BO for event energy flow at LEP1 and LEP2. 



Comment: "All variables are in good agreement with TWOGAM predictions". 
Note that the px range of jets may be of order of < Q 2 > for small Q 2 samples. 



29 



LEP2 - <Q 2 >=34 (GeV/cf 
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Figure 29: The comparison of data and the MC prediction for < Q 2 > = 34 GeV 2 (LEP2). 
a) Energy of tagged particle; b) The invariant mass; c) The number of jets; d) The jet 
angle with respect to the tagged particle; e) for jets in the same hemisphere as the 
tagged particle; f) p 3 ^ 1 for jets in the opposite hemisphere. Curves show the MC predictions: 
GVDM+QPM+RPC (GS2) (solid line), GVDM+QPM (dots), GVDM (dashes) (from EH). 



DELPHI LEP1 and LEP2 



■S 

-3 

2 25 
* 

-O is 

a, 

? » 



<(f>=13 (GeV/c 2 ) 2 




















--i m. . i . . . , 


i -3 -2 


-/ ! 


2 3 




T]=-ln(tan(®/2)) 

Figure 30: The event energy Bow as a function of the pseudorapidity. The notations as in 
Fig. || (from 



•OPAL 94 J2§ (LEP 1) 

Data were collected in 1990-92 for averaged Q 2 =5.9 and 14.7 GeV 2 . Early analysis 
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of the final hadronic state (the TWOGEN generator with the contribution based on 
QPM, VDM and on the FKP approach to the QCD contribution) was performed. For 
the estimation of the p? see Sec. 2. 1.2. 



•OPAL 97a @ (LEP 1) 

The hadronic final state was analysed in the measurement of the F% with one photon 
highly virtual (Q 2 between 6-30 GeV 2 - low Q 2 sample and between 60-400 GeV 2 - 
high Q 2 sample), the other being almost real. 

For generating the hadronic final state MC programs HERWIG, PYTHIA, and for 
comparison F2GEN were used with GRV and SaSID parton parametrizations. 

The energy (E t)OUt ) transverse to the plane defined by the beam axis and the tag 
direction, and other quantities for the low Q 2 sample are presented in Fig. 31. The 
discrepancy found for the E vis and E t)OUt distributions in this sample (Figs. |31] b and 
d, respectively) is absent in the high Q 2 events (not presented). 
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Figure 31: Comparison of data event quantities in the low-Q 2 sample with HERWIG, 
PYTHIA and F2GEN (PL) Monte Carlo simulations, a) the distribution of the visible 
invariant mass; b) the total visible energy of the event; c) the transverse momentum of the 
event in the tag plane; d) the energy out of the tag plane (from f24 



The hadronic energy E ttOU t distributions in the three XBj bins are shown in Fig. [32 . 



The failure of the models in the low Q region is most visible at low xbj- 

To establish the source of the discrepancy the energy flow per event in the low 
Q 2 sample was also studied as a function of the pseudorapidity (see Fig. |33[). The 
distribution of pseudorapidity r\ for the low Q 2 sample, corrected for the experimental 
effects, is shown in Fig. R3. 
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Figure 33: The hadronic energy Bow per event as a function of the pseudorapidity r\ for the 
data and various MC simulations, in various ranges of x v i s and Otag for the low Q 2 sample 
(from 



32 




Tl=-ln(tan(8'/2)) 

Figure 34: The measured energy flow in the low-Q 2 sample corrected for the detector 
inefficiencies, as a function of the pseudorapidity n, compared to the values generated by 
various MC models (from J2^/). 



Comment: None of the generators represents the final state accurately; E vis ,E t)OUt - 
distributions as well as the hadronic flow per event show a clear discrepancy. The 
failure of the models in the low Q 2 region is most marked at low XBj- The differences 
between MC models and data in the low Q 2 region are apparent, when the energy 
flow per event is plotted as a function of the pseudorapidity and the azimuthal angle 
(not shown). "Particular attention will need to be given to the angular distribution of 
partons in 7*7 system". 

The relation between p\ and Q 2 in the low Q 2 sample may indicate a need to take 
into account the structure of the virtual photon. 

•OPAL 97b (LEP 1) 

The measurement of F% and the modeling of the 7*7 fragmentation into hadrons at 
low Q 2 region (1.1 to 6.6 GeV 2 ) and very small Xgj - bins from 0.0025 to 0.2 is re- 
ported. The hadronic energy flow as a function of rj for three xsj regions is plotted in 
Fig.0 

Comment: Differences in the energy How distributions versus the pseudorapidity 
(Fig. |35[ ) and in the summed energy transverse to the tag plane are found for x vis < 0.05 
between the data and MC models (HERWIG, PYTHIA, F2GEN). 

•OPAL 97d H|] (LEP 1) 

An analysis of the hadronic final state was done, in which the discrepancies between 



the data and predictions, reported in OPAL 97a p4| , were examined from the point 



33 



of view of the number of produced jets. The data for Q 2 w 6 — 30 GeV 2 taken in the 
years 1994-95 were compared to the results of the HERWIG, PYTHIA and F2GEN 
generators. 
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Figure 35: The measured hadronic energy flow as a function of the pseudorapidity: a 
comparison of data with various MC predictions for three bins d two bins in 9 tag 

(from J||/j. 



The numbers of events for different number of jets, divided by the sum of all events, 



are presented in Fig. 36 
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Figure 36: The fraction of events with 0-3 jets. The points show the data with the back- 
ground subtracted, with statistical errors; histograms obtained with HERWIG, PYTHIA 
and F2GEN generators (from l3(ti ). 



The results of further studies of the energy flow versus the pseudorapidity 77 for 
events with different number of jets are shown in Fig. pTj. 
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1 jet events 2 jet events 

Figure 37: The hadronic energy Bow per event in bins of the pseudorapidity rj. Plots a), b), 
c) and d) show the average energy how per event for summed events, and for events with 0, 
1 and 2 jets, respectively. The samples are represented as in Fig. (from pfij/J, 



The number of events versus the energy transverse to the tag plane, obtained by 




E T [GeV] 
2 jet events 

Figure 38: The number of a) 1-jet and b) 2-jet events as a function of the jet transverse 
energy, E J T et . The HERWIG, PYTHIA and F2GEN events are normalized to the a) 1-jet 
and b) 2-jet data events (from J36j/). 
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Comment: "There is a marked difference between the data and the Monte Carlo sam- 
ples in the number of events with 2 jets' 7 . "While the F2GEN sample, generated using 
the pointlike approximation, overestimates the 2 jet rate by a factor 2.4, HERWIG and 
PYTHIA are too low by factors of 2.4 and 3.6 respectively." 

Disagreement between the data and each of generators for the hadronic energy flow 
versus the rapidity is seen for all types of events: 0, 1 and 2- jet. 
"All of the Monte Carlo samples model well the E tout distributions for events with 1 
jet, but PYTHIA underestimates the E t out of the events with 2 jets. " 



•OPAL 97e |J7fl (LEP 1) 

An analysis of OPAL data taken in 1993-95 was performed at < Q 2 > = 7.5 and 14.5 
GeV 2 . The possible improvement of the description of the transverse energy flow in 
DIS 7 experiment by introduction of the additional intrinsic k t smearing, influencing 
the angular distributions of the hadronic final state in the HERWIG, PYTHIA and 
ARIADNE generators was studied for the first time. 

In Fig. BP] the transverse energy out of the tag plane is presented together with the 
results of the simulation by PYTHIA and HERWIG generators. The energy flow as a 
function of rj was also studied for different x bins and Q 2 values (not shown). 
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Figure 39: The transverse energy out of the tag plane compared to PYTHIA (left) and 
HERWIG (right) simulations (from fflQ). 



Comment: "Both the E tjOU t and the energy flow per event are greatly improved with 
the inclusion of the power-like distributions of the intrinsic k t , with the exception of 
the peak in the energy flow at low x V i S and Q 2 = 14-7 GeV 2 , which still falls short of 
the data. " 

•TOPAZ 94 (TRISTAN) 

The jet production (one and two jets) has been studied in the deep inelastic e^y scat- 
tering (3.0 GeV 2 < Q 2 <30 GeV 2 ), see also Sec. 2. 1.2. Events with the transverse 
momentum of jets between 2 and 8 GeV and for the \rj\ < 0.7 were studied using the 
jet cone algorithm with R=l (see next section for the definition of the jets). 

The point-like and hadron-like configurations resulting in the different final-state 
topologies were studied. 
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Comment: The two high pr jet events are consistent with the point-like perturbative 
part. "In the one - jet sample an excess over the point like component is observed, 
which is direct experimental evidence for the existence of the hadron-like component in 
DISe-f". 
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2.3 Large processes in 77 and "fp collision 

The structure function F% considered in Sec. 2.1 is sensitive mainly to the combination 
of quark densities, and moreover, due to strong dependence on the charge, to the up- 
type quark distributions. Therefore it is of great importance to have an additional and 
very different from DIS 7 -experiments source of information. Hard processes involving 
resolved real photons provide such a source, with the cross section being a combination 
of the different parton densities convoluted with the cross sections for basic partonic 
subprocesses (see below). The role of these processes in determining e.g. a gluon 
density in the photon is unique indeed. 

As we have already mentioned, the jet production in 77 collision leads to the 
complementary information on the photon structure to that coming from the K7; 
moreover, the study of the hadronic final state has become recently a standard part 
of the structure function analyses performed at LEP (see discussion in Sees. 2.1 and 
2.2). 

In this section we focus on the large px jet production in the resolved 77 and •yp 
collisions, where photon may interact through its partons. The main goal here is to 
extract the individual parton density in the real photon 0. Beside the jet production 
also the inclusive one-particle production is sensitive to the partonic content of photon, 
but it depends on the additional fragmentation functions, and will not be discussed 
here 0. Also, other hard processes with resolved photon(s), like the Drell-Yan pair 
production may be used in pinning down the parton distributions in the photon. The 
lack of data reflects much lower rates for these processes. 

The production cross section of jets, which are the hadronic representations of the 
hard partons produced in the basic subprocess, does depend however on the applied 
jet definition. Although this introduces additional uncertainty in the description of 
data, relatively high rate for these events makes large pt jet production, apart from 
the F2 , the basic source of information on the partonic structure of 7. 

In the previous section we have discussed single tagged events with (positive) mass 
of the photons radiated by the tagged electron larger or much larger than 1 GeV (with 
exception of one experiment DELPHI 95). In this section the initial photons may 
be considered real , as the transverse momenta of the observed jets are chosen to be 
much larger than the median mass of the photon(s). 

At present, the 77 collisions arising in e + e~ colliders in no-tag Q conditions corre- 
spond to the 77 events where (both) real photons can be described by the Weizsacker - 
Williams energy spectra. In the OPAL experiment, the squared target photon masses 
were estimated to be < P 2 2 >=0.06 GeV 2 , in the ALEPH experiment around 0.23 
GeV 2 . In the 7^ scattering at ep collider HERA similar Weizsacker - Williams spec- 
trum describes the flux of photons coming from the electron (at HERA at present, i.e. 
in 1996/7, the positron) - here P 2 < 4 GeV 2 (but with the median 0.001 GeV 2 ) or 

10 Presently most of the data deal with the almost real photon; for the jet production in processes 
with virtual photon, see next section. 

n with the exception of the newest charged particle production measurements by HI group, where 
the gluon content has been derived. 

12 The name "untagged" is also used. In some cases antitagged or single - tagged events in the large 
Pt jet production are discussed in 77 and "fp collisions as well (see below). 
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below 0.01-0.02 GeV 2 , if the dedicated detector is used. 
2.3.1 General framework 

Resolved photon processes with large px jets can be characterized by the scale of 
hardness Q 2 which is usually provided by the px of the final jets, i.e. Q 2 ~ p\. (We 
introduce here the notation Q 2 in order to distinguish it from the DIS scale Q 2 , which 
is equal to the virtuality of the photon probe.) In order to resolve both photons the 
scale Q 2 should be much larger than P 2 - the virtuality of the most virtual photon in 
the process, and moreover much larger than Aq CD in order to apply the perturbative 
QCD. (Events where the masses fulfill the relation Aq CD C P 2 < Q 2 are discussed in 
the next section, where the concept and data on the structure of virtual photon will 
be introduced.) 

In this section initial photon(s) are considered to be real. Assuming the factoriza- 
tion between the hard subprocess and parton densities, the generic LO cross section 
for the two jet production in 7p collision (Fig. |40|a) is given by 

da*-** jet > x = J2 J J dx 7 dx p f ih (x 7: Q 2 )f j/P (x p , Q 2 )a»^ ^\ (19) 

where /i/ 7 (%, Q 2 ) describes the probability in the LL approximation of finding a par- 
ton of a type "i" in the photon. The x 7 (x p ) variable is by definition the part of the 
four momentum of the photon (proton) carried by its parton. Scale Q 2 is taken usu- 
ally to be equal to the so called factorization scale. Note that the precise definition 
of the factorization scale can be given in the QCD calculation going beyond the LO 
approach. 

For the two-jet production in 77 collision, 77 — > jet\ jeti X (Fig. ^0|b) a similar 
formula holds: 




Figure 40: Inclusive large px jet production in the resolved ■jp collision (a) and in the 
(double) resolved 77 collision (b). 

Note that, contrary to the DIS 7 experiments, in the resolved photon processes 
there is no one basic observable, analogous to the F%. Usually the pt or E T and 
pseudorapidity r\ distributions are compared with the QCD prediction using a specific 
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parton parametrization and the agreement justifies a posteriori the correctness of 
the applied parton parametrization and the QCD calculation. Only in few cases the 
parton distribution is extracted from data (see below). 

It is worth noticing the difference between the variables describing the produced 
parton and the corresponding quantities for the jet, representing the considered parton 
on the hadronic level. The unfolding procedure is needed in order to obtain the 
"true" partonic variables, e.g. transverse momentum of the parton (p t ) from transverse 
momentum or energy of the jet (j>t or Ex). 



2.3.2 Measurements of the resolved photon processes 

For the jet production processes where at the Born level (i.e. in the LO) two jets 
appear, the measurements of basic distributions in Et or the pseudorapidity r\ are 
performed for single or double jet (dijets) events as a rule. 

For comparison with QCD not only the above distributions are important but also 
the study of the structure of the jet (jet profile or transverse energy flow around the 
jet axis) and the energy of the underlying event, or the hadronic activity outside the 
jets where effects due to accompanying remnant jets may be seen. 

For more detailed study of properties of the resolved photon processes the sepa- 
ration of the direct photon events, where the photon participates directly in the hard 
subprocess, x 7 ~ 1, and resolved photon events is needed. Strictly speaking this sep- 
aration holds only in the LO approach. In the two jet events, x p and x 1 ( or x*' 2 ) 
distributions, sensitive to parton densities, can be reconstructed, see e.g. ZEUS 95b, 
97b. The angular distribution (e.g. in 77 CM system) da / d cos 9* , on the other hand, 
is not sensitive to these ingredients. So this measurement may help to verify the ex- 
pectation of different angular distributions which correspond to direct and resolved 
photon contributions (i.e. due to different subprocesses) , see ZEUS 96a for the first 
results. 

The theoretical predictions based on the LO |M| - or NLO(HO) QCD cal- 



culations are available for the inclusive one and two jet production for e + e and ep 



experiments |97j - ||110|| . The main MC generators used by the experimental groups 
for these kinds of processes are PYTHIA, HERWIG and PHOJET (with or without 
the multiple parton interaction). They allow to study different types of distributions 
of the final state hadrons. 

The following two jet definitions are used in the analyses reported below. The 
jet cone algorithm with the fixed value of the cone variable R, defined as R = 



(5(f)) 2 + (drj) 2 , with 5<p (5rj) describing the differences between the cone (jet) axis and 
the particle direction in the pseudorapidity and azimuthal angle, is used in PUCELL 
and EUCELL generators. The second one corresponds to the /^-cluster algorithm 
on which the KTCLUS generator is built up. Beside R also R sep is included in some 
analyses, which corresponds to the additional separation between partons - note the 
R sep = 2 R means no restriction. (The discussion of jet definitions can be found in 



13 



See the collection presented in Appendix. 
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The basic distributions, e.g. the jet Et distribution, are in general in agreement 
with the expectation for both single and double jet production. In the pseudorapidity 
dependence the discrepancy is observed both for the jet rates and for the transverse 
energy flow around the jet axis, especially for small Et and a small x 1 . This is taken as 
a hint that the multiple scattering may be important in the photon induced processes. 
Since the direct events should be free from such multiple interactions, the 77 and 7p 
collisions may offer a unique laboratory to study this problem. The strong dependence 
on the proper choice of the jet definition and the parameters like the cone size R, is 
observed recently in the jet production at HERA. There is a possible relation of this 
effect to the problem of describing the underlying events observed both in e + e~ and 
in ep collisions. 
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2.3.3 Jet production in 77 collision 



Here the results for the jet production in 77 collision, where one or two resolved 
photons may interact, will be presented. 

As in previous section we will denote the squared (positive) virtualities of two 
involved photons by P\ and P 2 2 - I n e + e~ machine parton momentum fractions - 
variables determined from the two final jets with the highest Et 

according to the formulae: 

x ± = Z jets (E±p z ) 

^hadrons(E ± p z ) 



DATA 



•DELPHI 94 HH (LEP 1) 

Data were collected in the period 1990-92 in antitagging conditions [P^ax ~ 0-12 
GeV 2 ). Jets defined according to the Lund cluster algorithm with pt greater than 
1.75 GeV were observed. 

Results for the pt of the jet distributions are shown in Fig. [|l] (p T has been intro- 
duced to separate the resolved photon contribution calculated using the FPK approach 
with DG and DO parametrizations). 

DELPHI 




2.5 3 3.5 
£ P T Jet(GeV/c) 

Figure 41: Distribution of jet transverse momentum compared with MC predictions 
(TWOGAM generator). Dark area - QPM only; dots - QPM + VDM; full line - QPM 
+ VDM + DG (pg, = 1.45 GeV/c); dashed line - QPM + VDM + DO (p° T = 1.22 GeV/c); 
(from E3I). 



•DELPHI 95 |3§] (LEP 1) 

(See also previous section). The values of Et of observed jets were larger than 1.5 
GeV, while the magnitude of the mass of virtual photon was equal to < >~ 0.06 
GeV 2 , i.e. the standard resolved (almost real) photon process, from the point of view 
of both photons. 

Comment: This is not typical DIS experiment since the photon probe is almost real. 
The analysis is not typical, either, for the standard large px jet study. 
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•DELPHI 97 (LEP 1 and 2) 

Based on the LEP1 and LEP2 data (1995 and 1996 runs), various distributions for the 
hadronic final state in no-tag events were measured (Fig. [42]) . 





Invariant mass, GeV/c 2 p/e<, GeV/c Pseudorapidity 

Figure 42: Distributions of: a), d) - invariant mass, b), e) - p^ ts , c) number of reconstructed 
jets, f) number of events versus pseudorapidity. Solid and dotted lines - VDM + QPM + RPC 
with GS2 and GRV parton densities, respectively. Dashed line - QCD - RPC contribution 
for GS2; hatched histogram - VDM part, double hatched - QPM part (from j^). 



Comment: "MC prediction (based on TWOGAM generator) gives perfect agreement at 
CM energy = 91 GeV and slightly exceeds data at 172 GeV ". 



•OPAL 97f (LEP 1.5) 

The inclusive one- and two-jet cross sections have been measured at ^fs~rZ =130 and 
136 GeV (based on the 1995 run). The anti-tagging condition was applied to one of 
the initial photons corresponding to maximum photon virtuality ~ 0.8 GeV 2 . 

The cone jet finding algorithm with R=l was used for the first time in photon - photon 
collisions at LEP. 

The transverse energy of jets Et is taken to be larger than 3 GeV and the pseu- 
dorapidity lies within \rj\ < 1. Analysis of the transverse energy flow around the jet 
direction for all jets and in two-jet events was performed. 

For the two-photon events, the distribution and the transverse energy flow 
around the jet, studied separately for the direct (x 1 > 0.8) and the double resolved 
(x 1 < 0.8) photons contributions, were obtained (see Fig. [|3]and Fig. [P], respectively). 
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Figure 43: The number of two-jet events as a function of x+ compared to PHOJET (solid 
line) and PYTHIA (dashed line) simulations. The hatched histogram is the direct contribu- 
tion to PYTHIA events (from Mj). 




Figure 44: The transverse energy Bow around the jets in two-jet events. The data are 
compared with the PHOJET (solid line) and PYTHIA (dashed line) simulations (from fjCj/). 



The one-jet and two-jet cross sections 



da 



and -fj^t were measured up to Ej, 



dE; 



jet 



jet_ 



16 GeV, extending the previous measurement at TRISTAN. Results are presented in 



Fig. |45 
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Figure 45: The inclusive one-jet (a) and two-jet (b) cross section as a function of for 

The solid line is the sum of 



the jets with \rji et \ < 1 compared to the NLO calculations \10t 
direct, single - resolved and double - resolved contributions shown separately (from p^/j. 
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The T) distributions for one and two jet events were studied as well (see Fig. [i"6|). 
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Figure 46: The inclusive one-jet (a) and two-jet (b) cross section as a function of E 3 ^ 1 
for the jets with > 3 GeV compared to the LO QCD calculations of PYTHIA 

and PHOJET generators (from J4?]/ ) . 
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Comment: "The data on rj distributions agree well with NLO QCD calculations based 
on GRV parametrization. The GRV-LO and SaS ID parametrizations describe the 
data equally well, the LAC 1, however, gives twice the observed value" . 



•OPAL 97g gl| (LEP 2) 

The dijet production in two photon collisions at e + e~ energy 161-172 GeV for the 
average squared virtuality P 2 = 0.06 GeV 2 was measured. The transverse energy of 
jets is taken to be Et > 3 GeV and the pseudorapidity lies within \r]\ < 2. The cone 
jet finding algorithm with R=l was used. 

The direct and resolved photon subprocesses (x^ < 0.8) were studied. The trans- 



verse energy flow around the jet axis is presented in Fig. £47] for the double resolved 
(a) and direct (b) photon processes. Larger hadronic activity for the resolved photon 
sample is observed as expected. 
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Figure 47: The transverse energy Bow around the jet axis in two - jet events; a) double 
resolved events, (b) direct events (from JS)/). 



The angular dependence plotted separately for the three basic resolved photon 
subprocesses and for the direct events is presented in Fig. |48. 
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Figure 48: The angular distribution of events separated into the direct and the double- 
resolved contributions, together with the QCD expectations (from p^/j. 



The cross section for the direct, single- resolved and double-resolved two-jet events 
versus Et compared to the NLO calculation with the GRV-HO parton parametrization 
is shown in Fig. f|9[ 
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Figure 49: The inclusive two-jet cross section as a function of E^* for jets with \r]^ et \ < 2, 



compared to the NLO calculations Solid line is the sum of the direct, single-resolved 

and double-resolved cross sections shown separately (from |H|/ ). 



The jet pseudorapidity distribution is presented in Fig. |50| together with the predic- 
tions from PYTHIA and PHOJET generators, with various parton parametrizations. 
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Figure 50: The inclusive two-jet cross sections as a function of \rf et \ for jets with E 3 ^ 1 > 3 
GeV (from g^j). 



Comment: "The Ex dependent two-jet cross section is in good agreement with NLO 
QCD calculation. The GRV-LO and SaS-lD parametrizations describe the inclusive 
two-jet cross section equally well. The LAC1 parametrization overestimates the inclu- 
sive two-jet cross section significantly" . 

•AMY 92 |@2l (TRISTAN) 

The measurement of the high px hadron production in the quasi-real 77 collision was 
performed (for the energy between 55 and 61.4 GeV). In the observed 3- and 4-jet 
events there are one or two spectator jets coming from the resolved photon(s). Two 
paramerizations of parton density in the photon were used: DG and (DO + VDM) 
(Fig. [H]) ( The p T has been introduced to separate the QPM+VDM contributions 
from the "MJET" one, which corresponds to the RPC contributions). 

500 
400 
300 

T3 

S 300 

IOC 



"0 2 * 6 B 

Pf- (GeV/c) 

Figure 51: The experimental p^ 1 distribution compared with the predictions of QPM (dot- 
ted histogram), QPM + VMD (dashed histogram) and QPM + VMD + MJET (solid his- 
togram) with pj, = 1.6 GeV (from p^/j. 
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Figure 56: Inclusive a) jet and b) two-jet cross sections as a function of px for the 
pseudorapidity \ r]j e t |< 0.7. TOPAZ Collaboration results from 



Comment: The polar angle distribution gives .."direct evidence of the presence of the 
resolved processes. The data exclude the parametrizations of LAC 3 and DO+ VMD, 
which predict a very large gluon content even at large x 7 ." 
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2.3.4 Jet production in the resolved 7p scattering 



Here we present the data for the jet production in 7p scattering taken at HERA 
collider at J~s^ ~ 300 GeV 0, where only one photon may be resolved [5 The 
photoproduction events correspond here to the limit of virtuality of the initial photon: 
P 2 ~ 4 GeV 2 (with median 0.001 GeV 2 ) or (with a special condition) ~ 0.01-0.02 
GeV 2 (see below). The partonic variables, both related to the initial and final states 
in the hard partonic subprocesses, are reconstructed from the corresponding quantities 
for final state hadronic jets. The relation between variables corresponding to these two 
levels depends on the order of the perturbative QCD calculation, as it was mentioned 
before. 

The parton momentum fraction x^ 1 (called also x v J s or x° bs ) is in the LO case equal 
to Xj. In practice it is reconstructed using the two jets with the highest transverse 
energy in the event, using the following relation: 

4* = — > ( 22 ) 



where also jet pseudorapidities, rf et = — ln(tan(6*/2)), and the energy of the photon 
E 1 (= yE e ) enter. (The positive pseudorapidity corresponds to the proton direction.) 
The scaled energy y of the initial photon (Eq. (^)) is measured from the transverse 
energy E^ and pseudorapidity of hadrons rj according to the formula 

» = 2jE^. (23) 

e h 

where the sum is over produced hadrons. 



In the following discussion the pseudorapidity of the jet in the laboratory system is 
denoted by rj, whereas rf means the corresponding variable in the 7£> CM system. For 
the HERA collider typically 77—77* ~ 2. The difference of the two jets' transverse energy 
Et\-Et2i AEt, and analogous difference for the pseudorapidities, A77 = 1 771 — 772 1 , are 
introduced in two-jet events. (For the analysis of the energy of underlying event the 
variable 5rj = r] ce u — 77 is used as in the 77 case.) In the analysis of two jet events there 
appear also the average pseudorapidity fj — (rji + r] 2 )/2 and the average transverse 
energy of jets, E T . 

Below, (f) denotes the azimuthal angle of the particle. In the cone jet algorithm, be- 
side R, also R se p is introduced both in the theoretical calculations and in experimental 
analysis of the jet production at HERA. 

Note also that the multiple interaction included in the experimental analysis is 
modeled as in the pp experiments (see e.g. HI 96a). 



14 The energy of the electron was at the beginning of running of the HERA collider equal to 26.7 
GeV, starting from 1994 it is 27.5 GeV, with the energy of the proton 820 GeV. 
15 The 77 events leading to the large pt jets are rare at HERA. 



51 



DATA 



•HI 92 g5[ (HERA) 

The evidence for the hard photoproduction of jet with E T > 10 GeV is reported (the 
jet cone algorithm with R < 1 was used in the analysis). 



•HI 93 i§ (HERA) 

First measurement of the inclusive jet cross section at the ep collider HERA (based 
on the 1992 data) is reported. The events with the scaled energy y of the initial pho- 
ton between 0.25 and 0.7 and the photon virtuality P 2 smaller than 0.01 GeV 2 were 
collected. The photoproduction of jet was studied for Et from 7 to 20 GeV and the 
pseudorapidity interval — 1 < r] < 1.5. The jet cone algorithm with R=l (PYTHIA 
5.6) was used. 

The transverse energy flow around the jet axis was studied and the discrepency 
was found in form of too large averaged E t on the forward side of the jet (see Fig. 
Among others, multiple parton interactions were mentioned as a possible explanation 
of this effect. 




Figure 57: The transverse energy Bow as a function of 5n (integrated over \5<j)\ < 1.0) (a, 
c) and as a function of 5(p (integrated over \Ar]\ < 1.0) (b, d). Figs, a) and b) correspond 
to -1.0 < i] < 0.5, b) and c) to 0.5 < rj < 1.5 (from BflJ). 



The inclusive jet cross sections da/dE T versus E T and da/dt] as a function of rj 
integrated over the corresponding range of the t] and Et, respectively, were measured 
and compared with the LO calculation using the following parton parametrizations: 
for the photon LAC2, LAC3 and GRV-LO, and GRV for the proton (see Fig. ||). 
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14 16 

E, [GeV] , 
Figure 58: Inclusive jet Et spectrum integrated over rj interval -1.0 < n < 1.5 (a) and 
inclusive n spectrum (b) for jets with Et > 7 GeV. LO QCD predictions generated by the 
PYTHIA generator using parametrizations LAC3 (dashed line), LAC2 (dashed-dotted line), 
GRV-LO (full line) and GRV-LO (without gluons, dotted line) (from g^). 



The shape of the do jdE^ is well described in the range of r\ between -1 and 1.5. It 
is a problem to describe the r\ distribution of jets. 

Comment: "In the (pseudorapidity) range 0.5 to 1.5 the data show larger average val- 
ues of the transverse energy flow outside the jet cone on the forward side of the jet 
than predicted by the MC. " 

LACS gives cross section higher by factor 3 than data for the Et distribution of the 
jets. "None of the models describe well the measured r\ dependence (for jets)." 

•HI 95 @7| (HERA) 

The photoproduction of 2-jet events (the 1993 data) was studied for the jet Et range 
from 7 to 20 GeV and the pseudorapidity interval < r] < 2.5 with |A^| < 1.2 (be- 
tween the most energetic jets) (see also HI 96a). The scaled energy y of the initial 
photon was between 0.25 and 0.7. Photon virtuality P 2 was smaller than 0.01 GeV 2 . 
The jet cone finding algorithm with R=l (and 0.7 for cross checks) was used (PYTHIA 
5.6 with the GRV-LO for the proton). 

For the first time the inclusive (LO) cross sections were derived for the parton level 
and the gluon density in the photon was measured. 

The transverse energy flow around the jet direction per event versus the rapidity 
distance from the jet direction was studied for 7 < Et < 8 GeV and < n < 1 and 
found to be asymmetric and different for the samples with x 7 > 0.4 and a; 7 < 0.4 (not 
shown) . 

The transverse energy flow versus the azimuthal angle around the jet direction and 
the transverse energy of the underlying events outside the jets (here named £]P edestal ^ 
another name - the underlying event) was studied (see Fig. |59|) with the conclusion 
"the multiple interaction gives an improved description" . 
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Figure 59: (a, c) The transverse energy How versus the azimuthal angle with respect to the 
jet direction in two rapidity bins, (b, d) Distributions of transverse energy measured outside 
of the jets. Histograms show PYTHIA simulations with (full line) and without (dashed line) 
multiple interactions (from p7|/). 



To achieve the goal which was here the extraction of the gluon distribution, the 
single parton cross section da/dp t integrated over the parton rapidity range, as well 
as the single parton da/drj were studied and compared with the LO parametrizations 
GRV and LAC 1 and LAC 3 (see Fig. ^0]) . 
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Figure 60: Single parton cross sections: (a) da/dpt integrated over the pseudorapidity range 
< 7] < 2.5, (b) da/dn for p t > 7 GeV. The solid line - the LO QCD calculation with 
GRV LO parametrization for partons in the proton and the photon. The dashed (dotted) 
line - the same for LAC1 (LAC3) parametrizations for the photon (GRV-LO for the proton) 
(from p^/j. 
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To extract the information on the subprocesses, the full two-jet kinematics was 
used and the distributions of the Arj, AE T , and of the x 7 and x v were studied. In 



Fig. 51 we present the distribution of the x 
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Figure 61: The two-jet event distribution of x 7 of the parton from the photon. The solid 
line - the contribution from the quark resolved photon processes; the dashed line - the direct 
photon contribution from the PYTHIA MC (from ^j). 



The (LO) gluon distribution in the photon was derived, at the average factorization 
scale < Q 2 >=< >=75 GeV 2 for 0.04< x 1 <1, see the table for numbers: 



<Q 2 > 


< x 7 > 


x 1 g(x 1 )/a 


[GeV 2 ] 




(stat. + syst.) 


75 


0.059 


1.92 ± 0.87 ± 1.68 




0.14 


1.19 ±0.34 ±0.59 




0.33 


0.26 ±0.24 ±0.33 




0.59 


-0.12 ±0.15 ±0.33 




0.93 


-0.08 ±0.61 ±0.30 



The results on x 1 g(x~ j )/a are presented in Fig. |62. 




Figure 62: The gluon distribution extracted from the resolved photon processes at 
< Q 2 >=< p 2 T >= 75 GeV 2 . For comparison GRV-LO (full line), LAC1 (dashed) and 
LAC3 (dotted) gluon parametrizations are shown (from p7(j). 



Comment: "The multiple interaction option gives an improved description" of the jet 
profiles and pedestal distributions ( still " deviation from data at large jet rapidities 2 
< 7] <2.5"). 

In the extracting the gluon density q 1 was taken as determined by two-photon exper- 
iments at LEP and KEK, in the form given by the GRV-LO parametrization. "A high 
gluon density at large parton momenta as suggested by the LAC 3 parametrizations is 
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clearly excluded. The strong rise of the LAC1 parametrization below x 1 < 0.08 is not 
supported. " 



•HI 96a gg| (HERA) 

The single jet production with Et > 7 GeV (and -1< 77 <2.5) was measured in ep 
collisions (data from 199?) with the scaled photon energy 0.25 < y < 0.7 and P 2 
below 0.01 GeV 2 . The cone algorithm with R=l was used. 

The properties of the hadronic final state and the distribution of the transverse 
energy are studied in detail (using PYTHIA, HERWIG and PHOJET generators with 
GRV-LO parton parametrizations for the proton and the photon). The integration 
over the 7p CM system pseudorapidity -2.5<r/* <1 leads to the total transverse event 



energy distribution shown in Fig. |63| a (here 0.3 < y < 0.7). 

The average transverse energy flow versus rj* for the total Ej- range between 25 
and 30 GeV was also measured (see Fig. |63|b). The shape of both distributions may 



indicate the need of the multiple interactions. 




Figure 63: a) The differential transverse energy cross section integrated over the pseudora- 
pidity (-2.5 < n* < 1). Histograms are the simulations with interactions of the beam rem- 
nants (full line - PHOJET, dashed - PYTHIA) and without them (dotted line - PYTHIA). b) 
The corrected transverse energy Bow versus rj* (n* > corresponds to the proton direction). 
The pseudorapidity range and histograms as in a) (from Ji^/). 



To get an insight into the details of the considered events the transverse energy 
flow outside of the two jets with the highest E T was studied as a function of x 7 for 
the |?y*| < 1 and Arj < 1.2. Results for the transverse energy density can be found in 
Fig. Bl 
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Figure 64: The corrected transverse energy density in the region \w*\ < 1 outside the jets, 
as a function of xt^ ts . The histograms are as in Fig. [6^ (from p^/j. 

The distribution of the transverse energy around the jet axis was also measured as 
a function of the 5<p (not shown). The jet width obtained in this analysis is similar to 
the corresponding quantity in the pp collision. 

Further results obtained for the jet cross section (here 0.25< y <0.7) are presented 
in Fig. |B~5|a, where da/dE? for the jet production in two rj regions is shown, and in 
Fig. [65|b, where the distribution da/dr] for the events with the transverse jet energy 
Et > 7, 11, 15 GeV is presented. The comparison with PHOJET and PYTHIA 
simulations, with and without multiple interactions, was done for both kinds of distri- 
butions. Note that the rapidity distribution is more sensitive to the photon structure 
functions, while the Et cross section to the matrix elements for the hard processes. 
Note also that "for Et bigger than 7 GeV previous measurements (HI 93) suffered 
from a defect and are superceded by this new measurement" . 

Comment: In addition to the primary hard scattering process, the interaction between 
the two beam remnants is included in the analysis. It gives "adequate descriptions 
of data" for the transverse energy versus pseudorapidity and the average energy flow 
obtained in this analysis. "For the first time the underlying event energy has been 
measured in jet events using direct and resolved photon probes. " 

The multiple interaction seems to improve also da/dET and da/dr] distributions for 
jets; within this approach the low Et and positive r\ range is still not properly described 
by the considered two LO parton parametrizations: LAC1 and GRV within PYTHIA 
program. PHOJET describes these data. 
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Figure 65: a) Cross section da/dE? for the jet production for E^ 1 > 7 GeV, in two rj 
ranges: -1 < rf et < 2 and -1 < rj jet < 1. The curves show the MC simulations with 
interactions of the beam remnants (full line - PHOJET, dashed - PYTHIA) and without 
them (dotted - PYTHIA). b) Cross section do/drf et vs rf et for different thresholds in Et- 
7, 11 and 15 GeV. The curves are as in a); additional dash-dotted curve - PYTHIA with 
LAC1 parametrization (from p^/j. 

•HI 96b (HERA) 

The double differential 2-jet cross section da / dx^ ts / d\og(E^ / Eq) was measured as a 
function of x 7 for different Ej- ranges above 8 GeV . Events (the 1994 data) correspond 
to: P 2 lower than 4 GeV 2 , the scaled energy y between 0.2 and 0.83 and the jet 
pseudorapidity -0.5< 77 <2.5, with At] < 1. The cone algorithm with R=0.7 was 
used. The PYTHIA generator with the multiple interaction using GRV-LO parton 
parametrizations for the proton and the photon was applied to describe the data. 

In Fig. |6^ the cross section is shown as a function of the x 1 for the different Q 2 
(= E^) bins (for the E\ distributions - see HI 97a). 
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Figure 66: The double differential two-jet cross section (see the text). The results from the 
PYTHIA simulation with multiple interactions and GRV-LO parametrizations are shown 
(from E|/j. 



For the first time the effective (LL) parton density of the photon 

9 , 



(24) 



with g 7 = E(^7 + ?t), was extracted, for 63 GeV 2 < p\ <1000 GeV 2 and 0.1< a; 7 <0.7 
(see HI 97a for the figures). 

Comment: "The effective parton distribution grows with the scale p\, although the 
increase appears slightly steeper than expected from the GRV-LO parametrization. " 



•HI 97a fl3H (HERA) 

A fixed cone algorithm with R=0.7 was used to describe jet events with P 2 < 4 GeV 2 
and y between 0.2 and 0.83. The pseudorapidity ranges < fj < 2 and Ar/ < 1 as well 
as AE T < 0.5E T were selected. 

The double differential dijet cross section data as a function of E\ for few ranges 



of x 1 are compared with the NLO QCD calculation KK fTUg and PYTHIA (GRV) 
simulation, see Fig. |67|. 
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Figure 67: The double differential cross section as a function of the square of the avera,. 
jet transverse energy and for the different xif ts (from 
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Figure 68: The effective parton density in the photon (from ffijj). 



The effective parton distribution in the photon was extracted for 0.2< x 1 <0.4 and 
0.4< x 7 <0.7. Its dependence on the Q 2 scale (= p\) is shown in Fig. |68. 



Comment: "Satisfactory overall description (of the double differential cross section 
for jet) except for x 7 > 0.6". 
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•HI 97b H5| (HERA) 

The new method of extracting the gluon density in the photon from the charged 
particles is introduced. Events with 0.3 < y < 0.7, P 2 < 0.01 GeV 2 and \r)\ < 1 were 
used. The result on the LO gluon density at Q 2 =< >=38 GeV 2 is presented in 



Fig. ^9| and compared with the jet data based on the 1993 runs (see also Fig. 61). 
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Figure 69: The gluon density (from j5|/). 
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•ZEUS 92 |g] (HERA) 

The evidence for the hard scattering in the photoproduction with E? > 10 GeV at 
HERA is reported (the jet cone algorithm with R = 1 was used). 



•ZEUS 94 Q53|] (HERA) 

The measurement was based on the 1992 data for the single and double jet photopro- 
duction for P 2 below 0.02 GeV 2 for tagged events, otherwise below 4 GeV 2 , and for 
the y between 0.2 and 0.7. The analysis of the direct and resolved photon processes 
was made using the HERWIG generator. The jet finding cone algorithm with R=l 
was used. 

The results for the E T distribution for single jets up to E T = 18 GeV, integrated 
over rapidity rj below 1.6, are presented in Fig. [?0|a. Fig. [70| b shows the da/drj data 
where the disagreement with MC prediction occurs for the positive rj. 
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Figure 70: Inclusive jet distributions for (a) transverse energy of jets, (b) pseudorapidity of 
jets (from J5^/J. 



Di-jet production has been studied by selecting events with two or more jets with 
Et > 5 GeV, for 77 smaller than 1.6. (Fig. fTl|). 
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rf* cos#" 
Figure 71: Kinematic distributions for events with two or more jets: (a) the jet pair invariant 
mass, (b) the transverse energy of jets, (c) the pseudorapidity, (d)cos8* of jet angles in jet-jet 
CM with respect to the proton momentum for events with > 16 GeV. The comparison 
with MC simulations is shown (from j5^/j. 



The x 7 and x proton distributions were studied as well for events with |A?7| <1.5, 
|A0| >120° and the invariant mass of two jets My larger than 16 GeV (not shown). 



•ZEUS 95a J5J (HERA) 

The 1993 data for the production of at least one jet with Ej- > 6 GeV are presented. 
Events correspond to P 2 below 4 GeV 2 , y between 0.2 and 0.85 and the jet pseudora- 
pidity between -1 and 2 (-3 and for rj*). PYTHIA 5.6 and HERWIG 5.7 generators 
(GRV and LAC1 parametrizations for the photon and MRSDo for the proton) were 
used with the cone algorithm for R=l. A wider than before range of Et (up to 41 
GeV ) and wider rj range were considered. 

The transverse energy flow around jet axis was studied. Results are presented in 
Fig. I72L where "there is some discrepancy for the forward-going jets in the 5r] > 1". 
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Figure 72: Transverse energy profiles as functions of 5n (top row) and 5(j) (bottom row) (see 
the text). Results from the PYTHIA simulation (with both resolved and direct processes) 
are shown (from |53/). 



The E T distributions integrated over two different pseudorapidity ranges and da/drj 
distribution integrated above three thresholds: 8, 11 and 17 GeV, are presented in 
Figs. 173, 173 for different parton parametrizations. 



Comment: "In the jet profiles, there is a significant excess of the transverse energy 
density in the data with respect to the MC expectations for jets in the region 1 < i] < 2. 
This excess is located outside of the jet in the forward direction Arj > 1." 
" Except for the region of very forward, low Et jets, these measurements are fully 
consistent with LO QCD predictions in new kinematical regime of the structure of the 
photon" . The result (for the da/dr) for the Et > 8 GeV and the range —l<r]<l) 
does not support the discrepancy of the da/dr] with respect to LO QCD calculations 
observed by the HI Collaboration (HI 93)". 
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Figure 73: The Et distributions for jets (see text). Results of the PYTHIA simulations 
with LAC1, ACFGP-HO, GS-HO and GRV-HO parametrizations of the parton densities in 
the photon are shown (from 



ZEUS 1993 



* 10 



10 



10 




GRV-HO 
Direct only (E, H > 17 GeV) 



ZEUS Data 
• E, w > 8 GeV 
■ if > 1 1 GeV 

V E," > 1 7 GeV 



"g 20 



\ 

b 10 



2 



. E,'" > 8 GeV 


, .i. I f fj 








f ZEUS Data 




/// 
/// 


LAC1 






ACFGP-HO(mc) 




GS-HO 






• • • GRV-HO 






Direct only 



2 



Figure 74: The jet pseudorapidity distributions: (a) integrated above three energy thresh- 
olds, Et > 8, 11 and 17 GeV (see text); the comparison with PYTHIA simulations using 
LAC1, GRV-HO parametrizations is shown; (b) only for Et > 8, and compared in addition 
with ACFGP-HO and GS-HO parametrizations (from |5i|/). 



•ZEUS 95b J5§ (HERA) 

The photoproduction of dijets, with at least two jets of E T larger than 6 GeV , is 
considered in the 1993 data. Events corresponding to the scaled energy y between 0.2 
and 0.8 and P 2 lower than 4 GeV 2 (for \Atj\ < 0.5) were grouped in the resolved and 
direct processes samples. The cone algorithm with R=l was used within the HERWIG 
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5.7 and PYTHIA 5.6 generators with the GRV-LO parametrization for the photon and 
the MRSD_ for the proton. 

The important x 1 distribution was studied. (The x 1 distribution has also been 
studied in ZEUS 96a, 96b, 97b.) The cut on the x 7 , equal to 0.75, was introduced 
later to enhance the resolved or the direct photon contributions, and a few distributions 
were studied separately for these samples. 



In Fig. [75] the transverse energy flow around the jet axis versus 5r) is shown, for the 
first time separately for the resolved photon and direct photon contributions ( ...with 
the failure to describe low x 1 data). 
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Figure 75: Transverse energy Bow around the jet axis versus Sn: (a) rr 7 > 0.75 (direct 
events); (b) x 7 < 0.75 (resolved events). Solid (dashed) line - PYTHIA (HERWIG) simu- 
lation (from pal). 



The da/drj was also measured for the direct and resolved photon events (not shown, 
see below for comments and new data in ZEUS 96a). 

Comment: "Both simulations fail to describe the transverse energy flow in the for- 
ward region (see also ZEUS 95a and HI 93). " 

"The LO QCD predictions (with DG, GRV and GS2 parton parametrizations) lie 
below the dijet cross section da/dfj data by factor 1.5-2". The importance of the NLO 
calculation is stressed. 



•ZEUS 96a |5§ (HERA) 

Analysis of the 1994 data for dijets (2 jets or more) for above 6 GeV and with 
the jet pair invariant mass above 23 GeV was performed. Events correspond to the 
range of y between 0.25 and 0.8 and P 2 below 4 GeV 2 . The cone algorithm with R=l 
was used within PYTHIA 5.7 and HERWIG 5.8 generators (with the MRS A parton 
parametrization for the proton and the GRV-LO for the photon). 

To obtain the scattering angle cos 9* distribution, sensitive to the parton dynamics 
and not parton densities as in analysis above, the cut not on A?7 (as in previous 
analysis) but on fj was introduced. 

Due to the 



The results for the x 7 , x p and 5r] distributions are shown in Fig. |76 



cut on r\ the absolute value of r\ is restricted to be below 1.8. Note that the applied 
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cut on the invariant mass suppresses events with low x. 
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Figure 76: Distributions in x p (a), x 7 (b) of the transverse energy How. The resolved (x 7 < 
0.75) and the direct (x 7 > 0.75) events as a function of Srj are presented separately (c and 
d, respectively). The MC results are also shown (from J5f|/). 

Very important results concerning the angular distribution due to various par- 
tonic subprocesses were obtained for the first time in the large px resolved photon 
processes. The angular distributions da/d cos 9* for the resolved and direct processes 
are presented together with the LO and NLO calculation based on the CTEQ3M 
parametrization for the proton and the GRV for the photon in Fig. [77]. The compari- 
son was also made with the HERWIG and PYTHIA simulations (not shown). 

Comment: The transverse energy flow is described properly, the " requirements of high 
mass and small boost remove the disagreement in the forward flow between data and 
the simulations which has been reported elsewhere in hard photoproduction at HERA. " 

The dijet angular dependence is well described by the LO and NLO QCD calcula- 
tions, and also by HERWIG and PYTHIA models. 
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Figure 77: The dijet angular distributions for the resolved and the direct contributions (see 
text) (from 



•ZEUS 96b |57[ (HERA) 

This is the extension of previous analyses (ZEUS 95, 96a) based also on the 1994 
data on the dijet production with P 2 lower than 4 GeV and y between 0.2 and 0.8. 
Results for the production of at least two jets with the pseudorapidity between -1.375 
and 1.875 and for E™ m = 6, 8, 11 and 15 GeV are presented, assuming |A^| < 0.5. 
In the data analysis different jet finding algorithms were applied: the cone algorithms 
EUCELL and PUCELL (both with R=l) and the k T - cluster algorithm KTCLUS. 

The resolved cross section was measured in the range 0.3 < x 1 < 0.75 and the 
direct one - for x 7 > 0.75. Analysis of the event distribution versus x 1 (see ZEUS 
97b for the figures), jet profiles in form of the transverse energy flow around the jet 
axis (presented in Fig. ^) and the do jdf] for the various jet definition and transverse 
energy thresholds (see ZEUS 97a,b for the new results) were performed. 
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Figure 78: Jet profiles < dEr/dSn > in different x 7 bins. Histograms based on the HERWIG 
generator with (solid line) and without (dashed line) multiparton interactions (from JSTj/j. 
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Comment: "inclusion of multiparton interactions improves the description (of the jet 
profiles) significantly." 



•ZEUS 96c [H] (HERA) 

The inclusive single jet cross section for the P 2 below 4 GeV 2 with y between 0.2 
and 0.8 and for three energy regions for W between 134 and 277 GeV was studied. 
The measurement of the transverse energy of the jets bigger than 14 GeV and the 
pseudorapidity range from -1 to 2 was performed. The HERWIG generator (using the 
jet cone algorithm with R=l) including the multiparton interaction was introduced 
in the analysis. The MRSA parton parametrization was used to describe the proton 
structure. 

Fig. [79] shows the rj distributions with the comparison to the LO (PYTHIA) and 
the NLO QCD predictions (KKS p3| ) (Fig. [fj a and b, respectively). 
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Figure 79: a) The differential cross section do /dif et integrated over the Ej? from four 
thresholds: E J T et > 14, 17, 21 and 25 GeV; PYTHIA results with the LO parametrizations 
MRSA for the proton and the GRV-HO and LAC1 for the photon are shown, (b) The 
same forE J T et > 17 GeV only; curves based on the NLO calculations (KKS \103j , using the 
GRV-HO and the GS-HO parton parametrizations) are displayed (from pBJ). 



The same is shown (for the NLO approach) in the form of the (data-theory) /theory 
plot for E T > 17 GeV in Fig. |(]. 

For events with Et > 14 GeV the rj distribution for different 7 — p CM energy 



ranges is plotted in Fig. 81 



Comment: "In the region r) > 1.5 the data show a flattening which is not described by 
the calculations. " 
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ZEUS 1994 Preliminary 
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Figure 80: The data as in Fig. [7|b with the prediction of the PYTHIA simulation (GRV- 
HO) and the KKS fcwSjl calculation (GRV-HO and GS-HO parton parametrizations) (from 
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Figure 81: The differential cross section do/dn^ 1 for three regions of the energy W. The 
PYTHIA results with the MRSA for the proton and with parton distributions in the photon: 
GRV-HO and LAC1 are shown (from ffil). 
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•ZEUS 97a [59| (HERA) 

The inclusive jet production da/dr] was studied with the iterative cone algorithm for 
events with P 2 < 4 GeV 2 and 0.2< y < 0.85. In the analysis R sep = R or 2R was 
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Figure 82: The differential cross section do /drf et integrated over Ej? from four thresholds: 
E J T et > 14, 17, 21 and 25 GeV; curves based on the NLO calculations KK^U§ using GRV- 
HO and GS-HO parton parametrizations for the photon and the CTEQ4M for the proton 
with a) R=l, b) R=0.7 are shown (from |5^/bj. 
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applied. The distributions of jets with E T lvn =lA, 17,21 and 25 GeV as a function of rj 
(-1< rj < 2) were measured. They are not properly described in the forward low 
region by the NLO calculation KK ||106|| if R—l is used, for R=0.7 the agreement is 
obtained (Figs. |82] a,b). 



The rj distributions in three regions of energy W are also in agreement with the 

1 (not 



NLO calculation for the R=0.7, what can be seen in Fig. The data for R 
shown) are not in agreement with a QCD calculation, as in Fig. 1ST 
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Figure 83: The differential cross section da/drf^ for three regions of the energy W. The 
PYTHIA results with the MRSA for the proton and with the parton distributions in the 
photon: GRV-HO and LAC1 are shown (from Mjb). 



Comment: "The measured cross sections for jets with R=0. 7 are well described by the 
(NLL QCD) calculations in the entire range ofrj." 

•ZEUS 97b J60[ (HERA) 

The dijet cross section based on the 1994 data for da/dfj for jets with Et > 6 GeV, 
1.375 < f] <1.875, and with the |A?y| < 0.5 was measured. The different Et thresholds 
and the jet definitions (as in ZEUS 96b) were implemented. The resolved and the 
direct photon events (0.3< x 7 < 0.75 and x 7 > 0.75, respectively) were separated in 
the analysis. The results together with the predictions from a NLO QCD calculation 
with the additional parameter describing the separation of jets R sep = R or 2R, as well 
as corresponding MC simulations (HERWIG 5.8 and PYTHIA 5.7 with or without the 
multiple interaction) are presented. 

The x 7 distribution obtained from the two-jet events indicates a need for the re- 
solved photon contribution (for the first time this kind of measurement was performed 
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in 1995, see ZEUS 95b and 96a,b). This analysis was based on the KTCLUS algo- 
rithm, for the result see Fig. 53. The small x 7 region is not properly described by the 
MC simulations both with and without the multiple interaction. 
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Figure 84: The corrected x° BS distribution, solid line - HERWIG with the multiple inter- 
action, dashed line - PYTHIA with the multiple interaction, dotted line - HERWIG without 
the multiple interaction (from 
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Figure 85: The uncorrected transverse energy How around the jet: solid line - HERWIG 
with the multiple interaction, dashed line - HERWIG without the multiple interaction (from 
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The transverse energy flow obtained using the KTCLUS algorithm is presented for 



different in Fig. |85|, with a similar discrepancy seen in the forward direction, as in 
previous measurements of the same quantity. 

The dijet cross section da/dfj obtained under condition \Arj\ < 0.5, for different 
Et thresholds and with different jet definitions was studied and the results are plotted 
in Fig. |S6|a. The predictions of the KK NLO QCD approach ||106|| are also presented. 

The same cross section, now with the KTCLUS jet definition, is plotted in Fig. |86lb. 
The predictions of the KK NLO QCD approach | |10(j|| with the different parton para- 
metrizations (GS and GRV) are compared with the data. 



ZEUS 1994 



3 3 
If 

I 2 
i 





WE. = 6 GeV i™ s >0.75 

/, at -v 1 

: EUCELL 
i PUCELL 
' KTCLUS 





.iineialerl l'nw 



0.5 



(b)E T "" = 8GeV j.™* 20.75 




|a:ne:pr.l eicel.ncJ error 



-1 



0.4 



(c)E r ~=110eV j^"" 5 > 0.75 



-Pi 



:i]iiaI|uI .iinvlaied error 



0.15 



0,05 



(d)E' "=l5GeV x, > 0.75 






ZEUS 1994 



WE,™ 


= 6GeV 


x™ S >0.75 






— GRVlKAK) 






GS (K&K) 






GRV (H&O) 






-j • ZEUS (KTCLUS) 








/A 






\ 








piinuip.; 


correlated erTor 



:o.4 - 



(c)E T ™=liGeV ^ > 0.75 




piriinp.J nerclaied error 



-1 


l 






i 




(flE I ~ = 8GeV0.3<x™ 


<0.75 


£ 






If 












e 


' It 


| t 




• 






principal correlated error 




-1 








1 \ 





4r 




(e)E, = 6 GeV 0.3< x, <0.75 




prainpd eterolated error 



F 

5 1 
D 



<».5 



£ M (g)E T ~= 11 GeV 0.3< i^"<B.15 

If 0.3 - 
B 

■° 0.2 - 

().] - 

o - 




pi; nuip.il airrdilL'J error 



S 0.1 

F 0.075 - 
B 

v 0.05 - 
0.025 - 
- 



(f)Ej™=8GeV 0.3< i,™ 


<0.75 




i 






principal correlated error 


-1 






1 n 



(h)E r mm =15GeV 0.3< ^ <0.75 



:'in:,i|i;il .widik'd uiTor 



Figure 86: a) Dijet cross sections obtained using the different jet algorithms. The curves 
show the results of the NLO calculation J J Off with R sep = 1 (solid line) and R sep = 2 (dashed) 
(GS?). b) Dijet cross sections obtained using the definite (KTCLUS) jet algorithm. The 



curves show the results of the NLO calculation |J0t^/ with R sep = 1 (solid line) and R 
2 (dashed) for different parton parametrizations (from pW). 



Comment: The discrepancy in the rj distributions for resolved photon contributions 
was found for events with E™ m > 6 GeV, as discussed in ZEUS 96b. 
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•ZEUS 97c 0] (HERA) 

The new 1995 data for the photoproduction of dijets (P 2 < 4 GeV 2 ) and y between 0.2 
and 0.85 are reported for the first time. The cone algorithm was used in the PYTHIA 
simulation; both jets have Et > 11 GeV and are in the rapidity range — 1 < r\ < 2. 
The preliminary data for da/dEx (symmetrized in 77; Et refers to the highest E T e ) 
are shown in Fig. [87] . 
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Figure 87: The dijet cross section da/dE? with one jet within the rapidity range 1 < r/i < 2 
and the other jet in three rapidity ranges, described in the figure (from /p^/J. 



The preliminary results for the pseudorapidity distributions are presented in Fig. 
for the direct and resolved photon samples. Below HERWIG and PYTHIA generators 
were used with the LL parton parametrizations for the photon. Similar analysis with 
the comparison to the KK | |106|| prediction is discussed in ZEUS 97d. 
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Figure 88: The rapidity distribution da/dr]2- The extra cut for the leading jet, Et > 14 
GeV. Three rji jet rapidity ranges were considered (as in Fig. 87). The comparison: PYTHIA 
and HERWIG predictions with the GRV-LO and the GS-LO parametrizations (from JS^/j. 
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•ZEUS 97d J62[ (HERA) 

The fey- cluster algorithm was used to study the dijet cross sections, 
d 3 cr / dE T dT]idr]2, for E T)1 > 14 GeV and E T>2 > 11 GeV and the pseudorapidity range 
-1< 771, 772 < 2 (see also ZEUS 97c). Note that y is between 0.2 and 0.85 as before, 
but P 2 < 1 GeV 2 . The Et dependence was also studied and compared with the KK 



106|| predictions (with the CTEQ4M and the GS96 parametrization for the proton 
and the photon, respectively). 

Below the results for the Et distribution in different regions of 771 2 are presented 



(Figs. |89]a, b). 
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The rapidity distribution da/dr) 2 integrated over Ej^ 
Fig. 00. Three r]i jet rapidity ranges were considered 



> 14 GeV is presented in 
(as in Fig. R3). Here the 



comparison was made with the NLO calculation by KK ||106|| with GRV-HO and GS96 
parton parametrization for the photon (CTEQ4M parametrization for the proton). 
The comparison with the HERWIG simulation was also made (not shown, see ZEUS 
97c). 

The high mass (M^ > 47 GeV) dijet cross section was also measured. Events with 
P 2 smaller than 4 GeV 2 were collected and compared with the same NLO calculation 
(KK ||106| ) as above. In this particular analysis the cone algorithm was used with R=l. 
Results are presented in Fig. |91l 
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Figure 90: The rapidity distribution do/dr\2- The extra cut for the leading jet, Et > 14 



GeV. Three rji jet rapidity ranges were considered (as in Fig. 87 ). The comparison with the 
NLO calculation, see text (from p2\] ). 
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Figure 91: The dijet distribution da/dcos9* for M%j > 47 GeV (upper plot) and dcr/dMij 
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3 Partonic content of the virtual photon 



The notion of partonic content of the virtual photon has appeared in high energy 
interactions soon after the related concept for the real one |], |13| . Not only the DIS 7 
formalism discussed before can be extended to the scattering on the virtual photon, 
but also the resolved virtual photon processes can be measured in this case. 

Note that for the resolved virtual photon processes the flux (convention dependent) 
of the virtual photon in the initial electron is introduced in the analysis |63|] (see also 
discussion below). The flux is usually taken (e.g. HI 97) in a form integrated over 
the relevant range of y, the scaled energy of the initial photon, and over the squared 
mass of photon, P 2 : 



Umax 



rn a x 



F l/e = J dyj p2 dP 2 f, /e (y,P 2 ) (25) 

V m in Tn.i n. 



with 



f~f/e{y,P ) = 7r-^[ w 2 — • 26 

2irP z y yP z 

For the photoproduction P^ in = m 2 y 2 /{\ — y), where m e is the electron mass. The 
assumption of the factorization of the cross section for the ep scattering and for the 
7*p process remains valid as long as p\ 3> P 2 . 



3.1 Theoretical framework 

The structure function of the virtual photon can be obtained in the Parton Model 
assuming the production of the qq pairs. In the mass parameter range 

m\ < P 2 < Q 2 (27) 

it has a form (to be compared with the Eq. @): 

a O 2 
F 2 7 (x, Q 2 , P 2 ) = N c Nj <Q A > -x{[x 2 + (1 - x) 2 } In ^ + 6x(l - x) - 2}, 

7T ± Ob 



where 



Nf 



<Q 4 >=^rY,Ql- (28) 



One can see clearly that the scale of the probe has to differ from the P 2 in order 
to test the structure of the virtual photon (the so called 'P 2 suppression' compared 
to the real photon case). The QCD evolution equations for the virtual photon are 
analogous to those for the real photon. Moreover in the case of the virtual photon 
there is a hope that the initial conditions are not needed, since for Q 2 ^> P 2 ^> Aq CD 



the nonperturbative effects should be absent (see Ref . fl3H ) . The virtual photon may 
play therefore a unique role in testing the QCD. (The existing parton parametrizations 
for the virtual photon are shortly discussed in the Appendix.) 
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The jet production in the processes involving one or two virtual photons can be 
studied in a similar way as for the real photons. The factorization between the emission 
of virtual photons and the jet production cross section for 7*7, 7*7* and 7*p scattering 
is usually assumed. Note that these cross sections involving virtual initial photons 
do depend on the convention used for the definition of the flux of virtual photons, 
since <r 7 * ~ l/T 7 */ e . At the same time the ee or ep cross sections are free from 
such ambiguity, as it cancels in the corresponding cross sections according to the 
(symbolical) relation: a e ~ T 7 . j e a 7 * . The notion of the resolved electron may happen 
to be very useful here (see the next section). 



3.2 DIS 7 * for virtual photons 

Measurements of the virtual photon structure function in the deep inelastic electron 
scattering are performed using double-tag events. Usually these events were selected in 
the kinematic region where one of the virtual photons (the probe) has, on the average, 
a large virtuality Q 2 and the other, the target, a small one, P 2 -C Q 2 . 
For the DIS 7 * measurements, the quantity 

x ms = Q 2 /{Q 2 + P 2 + W 2 ms ) 

needs to be converted to the true xsj = Q 2 /^Vl variable 0. We start with the old 
data from PLUTO, the only measurements of "DIS" type for the virtual photon so 
far. 

DATA 



•PLUTO 84 []64[ (PETRA) 

The double - tag events were measured where one of the virtual photons (the probe) 
had, on the average, virtuality < Q 2 >=5 GeV 2 and the other (the target) < P 2 >=0.35 
GeV 2 . (The energy of the beam was here 17.3 GeV.) The experiment was sensitive to 
the following combination of the virtual photon structure functions: 

F eff = F 2 + (3/2)F L . 



F e ff was extracted for XBj range between 0.05 and 0.6, see Fig. |92|. In Fig. £)3| the 
quantity Q z 'o " 77 / 'An" 2 ''a 2 : , averaged over both XBj and Q 2 , is shown as a function of the 
measured P 2 (0.2 - 0.8 GeV 2 ). 

16 At finite P 2 a modified variable Xgj, which extends over the whole range between and 1 may 
be introduced. For small ratio P 2 /Q 2 <C 1 they coincide. 
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Figure 92: The PLUTO Collaboration data for the effective structure function for the virtual 
photon. For the fixed averaged Q 2 = 5 GeV 2 and P 2 = 0.35 GeV 2 , a dependence on XBj is 
shown (from p4j). 
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Fi gure 93: The PLUTO Collaboration results for ^1^2 O77 f° r virtual photon as a function 
of its virtuality P 2 (averaged over xbj and Q 2 ranges) (from Jfej/j. 



3.3 Jet production in resolved virtual photon(s) processes in 
77 and 7£> collisions 

As for the real photon, the large pt jets may resolve the virtual photon(s). Provided 
the corresponding mass relation Q 2 ~ p\ 3> P X 2 (P 2 2 ) ^ ^-qcd occurs, one may use the 
QCD improved parton model as in the case of the real photon (Eqs. ( |I9"| , pT|)). 

Recently the LO and the NLO QCD calculations have appeared for the jet pro- 
duction by virtual photon(s) in the considered processes 77 and jp. 

The transition region between the interaction of an almost real photon and of a 
virtual photon with the proton is studied by HI and ZEUS collaborations in the ep 
collision at HERA. In such analysis the MC generators used to describe the photopro- 
duction and the DIS events as well as the rapidity gap events at HERA are used. As 
far as the flux of virtual photons is concerned, it is integrated over the corresponding 
range of the virtuality (see Eqs. (p6| ^271)). 
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Some kinematical variables are defined in the 7*p CM system, and are denoted 
below by a star, e.g. E^. Note that we use the notation P 2 for the squared virtuality 
of the photon although in the context of the DIS events (on the proton) at HERA it 
plays the role of the Q 2 . 

DATA 



•TOPAZ 94 Q2|] (TRISTAN) 

The jet production TOPAZ 94 data (see discussion in Sec. 2. 2) are related to the 
resolved virtual 7(P 2 )+real ^{P 2 ) process. The jet production (one and two jets) has 
been studied with 3.0 GeV 2 < P 2 <30 GeV 2 , and p T between 2 and 8 GeV (corre- 
sponding to 4< Q 2 <64 GeV 2 ). To what extent one can describe these events using a 
partonic language for the virtual photon is not clear since in some cases one is probing 
the region Q 2 /P 2 smaller than 1. (See also Sec. 2. 2.) 



•HI 97 |>§ (HERA) 

The single jet cross section for the events with > 4-5 GeV and with 0.3 < y <0.8 is 
studied in the transition between the photoproduction and the standard DIS P regime 
at HERA. The data from the years 1994 and 1995, in three ranges of the squared mass 
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Figure 94: The inclusive da ep /dE^ jet cross section as a function of the transverse energy 
Eji for various initial photon virtuality P 2 ranges and for the —2.5 < rj* < —0.5. The 
HERWIG prediction is denoted by the solid line, the dashed line corresponds to the direct 
contribution to this model (from J§^/ ). 
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of the virtual photon P 2 <1(T 2 GeV 2 , 0.65 < P 2 <20 and 9 < P 2 <49 GeV 2 were 
collected. The jet pseudorapidity was studied in the range —2.5 < rj* < —0.5, using 
the k T - clustering algorithm and the PHOJET 1.03, LEPTO 6.5 and ARIADNE 4.08, 
RAPGAP and HERWIG 5.9 generators (with the JETSET used for the hadroniza- 
tion). The GRV 94-HO parton parametrization for the proton and the Drees-Godbole, 
GRV-HO, SaS-2D parametrizations of the virtual photon were used. 

The measured da ep /dE^ as a function of the transverse energy of the jet for various 
P 2 ranges, integrated over y between 0.3 and 0.6, is presented in Fig. |93|. It was found 
to be in agreement with the HERWIG (DG) model. In Fig. [5^ the corresponding 
data for the rapidity distribution for jets with > 5 GeV are shown for the various 
virtuality ranges. 
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Figure 95: The inclusive da ep / dn* jet cross section as a function of the rapidity rf for various 
initial photon virtuality P 2 ranges for the > 5 GeV. The HERWIG (DG) prediction is 
denoted by the solid line, the dashed line denotes the direct contribution (from XBW). 



The results together with the HERWIG and RAPGAP predictions based on the 



DG and SaS-2D parton parametrizations in the virtual photon are shown in Fig. [96. 
To study the dependence of the virtuality of the photon the cross section <7 7 . p is 
introduced, 



a. 



•y*p->jet+X 



0~ep^jet+X 

~F~i 

1 7/e 



(29) 



although it is not certain that the above factorization really holds for the whole range 
of kinematical variables. 
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Figure 96: The inclusive j*p jet cross section as a function of P 2 (= Q 2 ) for various ranges 
of transverse jet energy E* T 2 = (Q 2 ) for the -2.5 < if < -0.5. The HERWIG (DG) 
prediction is denoted by the dashed line, RAPGAP (DG) - the dotted line, RAPGAP (SaS- 
2D) - the solid line and the dot-dashed line corresponds to the HERWIG with the GRV-HO 
parametrization as for the real photon ("no P 2 suppression") (from j^qj). 



Comment: The HERWIG (DG) model gives " a good description of the data except 
of jets in the lowest El range when 9< Q 2 < 49 GeV 2 ". 

•ZEUS 95c |§je|] (HERA) 

The measurement of the direct and resolved photoproduction at HERA with the vir- 
tual and the quasi - real photons was performed. In the 1994 run two samples of 
events were collected: with the photons of virtualities 0.1 GeV 2 < P 2 < 0.55 GeV 2 
and with the quasi - real photons (P 2 < 0.02 GeV 2 ). For each dijet event the fraction 
of the photon momentum x° bs , manifest in the two highest Ej> jets (Et > 4 GeV), 
was calculated. The events associated with the direct photon process (high x° bs ) and 
with the resolved photon processes (low x° bs ) were found both in the virtual and quasi 
- real photon samples. The ratio N res /N^i r was calculated as a function of P 2 and it 
seems to decrease with an increasing photon virtuality, see ZEUS 97c for new results. 

•ZEUS 97c |67|] (HERA) 

A study of the ZEUS dijet production data taken during 1995 for the transition region 
between the photoproduction and the DIS P is reported. In particular two samples were 
studied: 0.1 < P 2 < 0.7 GeV 2 and P 2 < 1.0 GeV 2 . The range of the scaled energy 
was 0.2< y <0.55. 
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The jets with E T > 6.5 GeV and -1.125< r] <1.875 were used. The HER- 
WIG 5.9 generator with or without the multiple interaction was applied (the MRSA 
parametrization was used for the proton and the GRV parametrization for the pho- 
ton). The events corresponding to the direct processes, with x 7 > 0.75, and to the 
resolved ones for x 7 < 0.75, were studied. 

The results for the x 1 distribution for the different P 2 ranges are presented in Fig. 



9% and the ratio o 'resolved/ direct versus P is plotted in Fig. |98 
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Figure 98: The ZEUS Collaboration results for the ratio of the events associated with the 
resolved photon to the events due to the direct photon processes as a function of the squared 
virtuality of the initial photon (from p7j). 
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4 Related topics 



4.1 On leptonic structure functions of the photon 

The lepton pair production by two photons in e + e~ scattering, being one of the basic 
QED processes^ is also a test for experimental methods applied to more complicated 
two-photon reactions involving hadrons. This applies especially to the process e + e~ — > 
hadrons where the hadronic structure functions of the photon are measured and the 
unfolding and tagging methods have to be tested (Sees. 2.1 and 2.2). 
The lepton pair production in single tagged events 

e + e" -> e + e~l + r (30) 

can be described by introducing the QED structure functions for photon: F^^ ED \ 
p^(QED) ^ p'yjQED) ^ ^ g ^ |68| ; |69|j). Unlike in the hadronic case {i.e. with the 
hadronic final state), these photon structure functions can be reliably calculated in 
QED. Moreover all the particles in the final state can be directly observed. 



Although the final state e + e — >e + e e + e has been also measured [72], for technical 
reasons only muonic structure function, denoted below as F^ ED > ( could be extracted 
(Eq. (|30| ) with I — p). 

pWJ ED ) s t rU cture function is much harder to measure, because its contribution 
is weighted by the small factor y 2 (as for hadronic F2, see Eq. (||)). However, this 
longitudinal structure function is not the only structure function that contains addi- 
tional information. It has been shown that there are azimuthal correlations in the final 
state particles from two-photon collisions which are sensitive to additional structure 
functions 0, [71]]. 

If, instead of measuring the cross section d °^ g " dy (Eqs. @, ©), one measures in 
two-body 7*7 collision also one final state particle a, additional structure functions 
FJ and Fg appear jTTJ: 



d(j(e7 -> eaX) 2na 1 + (1 - y) ~ ~ 

dxBjdydU a /A7r Q 2 x B jy 

-p{y)Fl cos a + \e{y)Fl cos 20J . (31) 

Here Q a describes the direction of particle a in the 7*7 rest frame, and cj) a is its 
azimuthal angle around the 7*7 axis, relative to the electron (tag) plane. The functions 
e(y) and p(y) are very close to 1. The standard functions Fj, and F^ are obtained from 
the corresponding F^ by integration over the solid angle Q a . Note that the formula 
(31) holds for two leptons or two partons produced in the final state. 

The function Fg ED ' is, in the LL approximation and zero muon mass limit, equal 
to F2^ E , although it involves quite different photon helicity structures. Thus ex- 
tracting jp^ E can give us indirectly information on F^ EC>S} (in LLA). 
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The earliest experiments on e + e — >leptons are summed up in |6SL 69 
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DATA (Early experiments) 



•CELLO 83 jgi (PETRA) 

Here both ee and uu pairs were observed. At < Q 2 >=9.5 GeV 2 the "muonic" struc- 
ture function F% was measured Q and was found to be in good agreement with 
QED calculations. 

•PEP-9/TPC 84 |73| (PEP) 

pi{QED) measured at < Q 2 >=0.3 GeV 2 . Results are shown in Fig. 105. The 
quantity 2x Bj F? {QED) was also extracted in this experiment, from the y region where 
its contribution was comparable to that of F2 ' . 

Comment: The importance of fi + n~ bremsstrahlung background subtraction has been 
shown in extracting F^ QED \ Reasonable agreement with QED predictions was ob- 
tained, especially for xsj < 0.25. 

•PLUTO 85 Q74| (PETRA) 

p^(QED) measured at < Q 2 >=5.5 and 40 GeV 2 in the full x B j range 18 . Results 
are in agreement with QED calculations. 

DATA (Recent results) 

In the last few years new measurements of the leptonic structure function p^ ED ^ 
have been performed at LEP by the four collaborations; the additional structure func- 
tions Fj^2 ED ^ have also been measured by ALEPH, L3 and OPAL groups. 

Since the invariant mass of the [i + pair, and hence XBj, can be determined 
very accurately, the measurements of F^ ED > are only statistically limited. This is in 
contrast to the hadronic final states, where the extraction of x B j introduces significant 
uncertainties. Thus the investigation of QED structure functions is no longer treated 
only as a test of QED, but rather as a clean experimental procedure meant for testing 
and refining the experimental procedures to be used in much more complex case of 
hadronic final states . 

Because of the precision of LEP data it is also possible to study the effect of the 
(small) virtuality P 2 of the quasi-real (target) photon. 

•ALEPH 97d J75[ (LEP 1) 

F J(QED) ^ been meagured ( data from lg94 ) for Q 6< Q2 <Q 3 Ge y2 Q2 > = 2 .8 

GeV 2 ) and 3.0< Q 2 <60.0 GeV 2 (< Q 2 >= 14.6 GeV 2 ) (Fig. M). 



18 Results for all F^^ EDS> measurements discussed here will be shown together in Fig. 107. 
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Results for 1994 Data Results for 1994 Data 
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Figure 99: Structure functions FrJ /a measured in the ALEPH experiment |7|/ at 
< Q 2 >=2.79 GeV 2 , < P 2 >=0.153 GeV 2 (a) and < Q 2 >=14.65 GeV 2 , < P 2 >=0.225 
GeV 2 (b). (from ffijjj). 



In this experiment also the azimuthal angle distributions have been measured and 
functions Fl' QED) and F^ QED) extracted (Fig. [100]). 
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Figure 100: The structure functions F fQ ED )/ F ^ ED ) ( a ) and i_ F 7(QED) /f ^qed) (b) 
measured in the ALEPH experiment for < Q 2 >=8.8 GeV 2 [75[/; the solid lines are the QED 
expectations, (from ffBj). 



Comment: Both F^ Er)S} an d p^ ED ) a g ree verv we n with QED predictions. Poorer 
quality of F^ ED ^ data can be helped with better statistics. 



•DELPHI 96a j$Z§ (LEP 1) 

p'jiQED) been measured at < Q 2 >=12 GeV 2 , as a test for the unfolding and 
tagging methods in extraction of the hadronic in DELPHI experimental environ- 
ment (see Sec. 2.1.2). The effect of non - zero target virtuality has been studied (see 
Fig. |101| ). A satisfactory fit to the measured F2^ ED ^ is obtained for the fixed value 
of P 2 =0.04 GeV 2 . 
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Figure 101: Unfolded F%^ E j a compared with QED predictions obtained for different 
masses of the target photon: zero mass (upper curve), < P 2 > = 0.04 GeV 2 (middle curve), 
and average value < P 2 > = 0.13 GeV 2 (lower curve) (from ^j). 



•L3 95 Q76fl (LEP 1) 

The data collected by L3 detector in the years 1991-93 were used to extract the F^ ED ^ 
structure function for 1.4< Q 2 <7.6 GeV 2 (< Q 2 >=1.7 GeV 2 ) (see Fig. gO^). In ad- 



dition, the angular distribution in the azimuthal angle a was measured and some 
information on previously unmeasured Fj^ ED ' and F^ ED > was obtained. 
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Figure 102: F^®**^ /a measured at < Q 2 >=1.7 GeV 2 . The effect of nonzero target 
photon mass is also shown (from |7q/ j. 



•L3 97 []77[| (LEP 1) 

The L3 Collaboration has performed a new measurement of QED photon structure 
functions: F2^ ED \ F^ E and F^ ED ^ in the e + e~ — > e + e~fi + {i~ process f32"f . 

The effect of the photon target virtuality has been studied and is clearly seen in 
all the three structure functions (see Figs. |103| and |104|) . 
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Figure 103: L3 data on F 2 7W ^ j (x, Q 2 , P 2 )/a compared to QED predictions at P 2 =0 (solid 
line) and < P 2 > fit = 0.033 GeV 2 (dotted line) (from frfij). 
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Figure 104: L3 data on structure functions F^ QED ^ /a and F]^ QED ^ /a measured at 1.4 
< Q 2 < 7.6 GeV 2 . The curves are the QED predictions for P 2 =0 (solid) and < P 2 >fu= 
0.033 GeV 2 (dashed), respectively (from |7"7|/). 



•OPAL 93 Rf§ (LEP 1) 

The QED structure function F^ iQED) for < Q 2 >=8.0 GeV 2 was extracted from single- 
tag events at CM energy ~ Mz- In Fig. |105| the data and QED expectations are 



presented together with the earlier measurements at CELLO 83 [72| and PEP 



9/TPC 84 [0]. 
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Figure 105: F^ iQED) /a measured at CELLO frSjj (open squares), PEP-9 (7|/ (open circles) 
and OPAL fl7q/ (full circles). The solid lines are the QED expectations (identical for CELLO 
and OPAL within the systematic errors) (from Wftl). 



•OPAL 97h []79[] (LEP 1) 

Here the extraction of F^ QED) was performed for 0.85< Q 2 <31 GeV 2 (< Q 2 > =5.2 
GeV 2 ). Fj(Q ED ' > was not extracted due to the partial integration of the cross section. 

The measured value of ^eF^ ED ^ / F^ ED ^ i s significantly different from zero and 
its variation with xbj is consistent with QED (see Fig. |106|) . 
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Figure 106: The values of ^F^ ED ^ / F^{QE D ) bt ameo ] f rom the azimuthal angle distri- 
butions (corrected for the effects of the detector). The solid line is the QED prediction for 
Q 2 =5.2 GeV 2 and e=l (from Mj). 



A compilation of results for the QED structure function F2 ' can be found 
in Fig. 107| (from P2"|). Data from old and new experiments, compared to QED 



calculations of F2^ QED \x, Q 2 ,P 2 = 0), are presented as a function of xbj for Q 2 from 
0.1-40 GeV 2 . 
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Figure 107: Summary of existing F% /a data for broad Q 2 range shown with QED 
predictions for P 2 =0 (from JS^/ j. 



4.2 On the structure function of the electron 

In the DIS experiments in e + e~ collisions the inclusive hadron production can be 
ascribed not to the photon - target but rather to the parent electron or positron 
- target. This in some cases may be more straightforward than pinning down the 
structure of the virtual photon, as was mentioned before. 

This topic is discussed in e.g. Ref. fjjf, @- See also the results in [jSOfl , where 
the structure function of the electron (in general - lepton) is related not only to the 
structure function of photon but also to the electroweak gauge bosons W and Z. The 
structure of weak bosons which appears in this approach was introduced and discussed 
in Ref. gl| . 
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5 Summary and outlook 



The status of recent measurements (~ 1990 and later) of the "structure" of unpolarized 
real photon in "DIS 7 " experiments as well as in large pt jet production processes 
involving resolved photon(s) is presented. A qualitative change has appeared in both 
types of measurements. Final results based on few years' runs at LEP 1 are being 
published with the higher statistics and improved unfolding methods. The new data 
from LEP 1.5 and 2 appeared during last year. On the other hand impressive progress 
has been obtained in pinning down the individual parton contributions in large px 
resolved photon processes both in e + e~ and in ep collisions. 

The existing discrepancies in describing the final hadronic states in the DIS 7 exper- 
iments as well as in the resolved photon processes in 77 and 7^ collisions seem to have 
the common origin. The need of additional p? in the distribution of produced hadrons 
and jets may suggest an extra interaction involving the constituents of photon(s). One 
of the possible explanations could be a multiple interaction described in a similar way 
as introduced in the pp processes. 

As far as the data related to the "structure" of the virtual photon are concerned, 
they have just appeared from the hard photon-proton and photon-photon collisions 
with resolved photon(s). The corresponding new DIS 7 * measurement of the structure 
functions of virtual photon will come probably only with the new generation of accel- 
erators. The interesting extension of the idea to other gauge boson structure functions 
and the related concept of the structure function of the electron might also be tested 
there. 

The new measurements of the partonic content of the photon are accompanied by 
the impressive progress made in the NLO QCD calculations for the resolved real and 
for virtual photon processes. 

The future high-energy linear e + e~ colliders LC as well as related ej and 77 collid- 
ers, based on the backward Compton scattering on the laser light, will offer a unique 
opportunity to measure the structure of photon in a new kinematical regime. More- 
over, at these colliders measurements of structure functions for photon with a definite 
polarization should become feasible with a good accuracy (see e.g.[p2|). 
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6 Appendix 



6.1 Parton parametrizations for the real photon 
Duke - Owens (DO) || 



A leading logarithmic parametrization of the parton distributions in an asymptotic 
form. Quarks with equal charges have the same distribution functions: /„ = f c , 
h = f s (JV/ = 4). 
Drees - Grassie (DG) @ 

A parametrization for a full solution of the leading order evolution equations. The 
input parton distributions with free parameters assumed at Qq — 1 GeV 2 and fitted 
to the only data on F% existing at that time, at Q 2 = 5.9 GeV 2 , from PLUTO. 
Field - Kapusta - Poggioli (FKP) [0 

In this approach F% is divided into the hadronic part (F^ 1 ) and the point-like one 
(F 2 PL ). The F^ L arises from the basic 7* — > qq coupling and higher order QCD 
corrections, if the final jet px is greater than p T . If Pt is smaller than p T then the qq 
pair creates a bound state. F^ AD in this non-perturbative case is taken from the VMD 
model. The perturbative _F 2 PL is calculated using the first order splitting functions and 
the one-loop as- 

Levy - Abramowicz - Charchula (LAC) |85| 

A parametrization for a full solution of the leading order evolution equation fitted 
to all available in 1991 measurements of F 2 for Q 2 > Qq. Three sets are provided 
with the different choices of an input scale Qq, and the x — > behaviour of a gluon 
distribution G(x), namely: 
• LAC1: Q 2 = 4 GeV 2 

Ql = 4 GeV 2 , xG(x) -> const. 
Q 2 = 1 GeV 2 . 



LAC2 
LAC3 



Gliick - Reya - Vogt (GRV) gB| 

The LO and NLO parametrizations of the parton distributions generated dynamically 
from the valence-like VMD input. The low initial scale Q^ = 0.3 GeV 2 is universal 
for the proton, the pion and the photon structure functions. The DIS 7 scheme is 
introduced to avoid the large-x^- instability problems. The one free parameter, which 
is a VMD input normalization, is fixed by the data. 



Gordon - Storrow (GS) |7 



The LO and NLO parametrizations. The input structure function at scale Qq = 5.3 



GeV |87|a and Qq = 3 GeV |87|b in the LO analysis is chosen as a sum of a hadronic 



part from the VMD model and of a point-like part based on the Parton Model. Free 
parameters (also light quarks masses) are fitted to the data for Q 2 > Q 2 . The NLO 
distributions in the MS are obtained by matching of the F 2 in the LO and the NLO 
approaches at the Qq scale. 

Aurenche - Chiapetta - Fontannaz - Guillet - Pilon (ACFGP) |S3]a 

A solution of the NLO evolution equation with the boundary condition taken at Qq 

= 0.25 GeV 2 . The input parton distributions was obtained from the VMD model at 

Q 2 = 2 GeV 2 and evaluated down to Qq. 

Aurenche - Guillet - Fontannaz flS8| b 
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The NLO parton distributions obtained with the input distributions (shown to be 
scheme-dependent) at Qq = 0.5 GeV 2 . The input distributions are based on the VMD 
model modified to agree with the MS scheme used in this analysis. 
Watanabe - Hagiwara - Tanaka - Izubuchi (WHIT) 

A set of six the LO parametrizations obtained by fitting the input distributions to all 
available data for 4 GeV 2 < Q 2 < 100 GeV 2 . The parametrizations WHIT1-WHIT6 
are based on different input gluon distributions. A massive charm contribution is 
calculated from the quark parton model and for Q 2 > 100 GeV 2 from the massive 
quark evolution equations. 
Schuler - Sjostrand (SaS) |J0 



Four sets of the LO parametrizations. The non-perturbative input distributions at 
Qo = 0.6 GeV in the SaSID, SaSIM and Q = 2 GeV in the SaS2D, SaS2M sets 
are based on the VMD model (their normalization is fixed and the x - dependence is 
obtained from the fits to the data). The fully calculable point-like contribution to the 
F% is expressed as an integral of the "state" distributions over the virtuality k 2 of the 
7* — > qq state. 

The non leading term C 1 is included into the F2, leading to the MS distributions 
(SaSIM, SaS2M). C 7 = gives distributions in the DIS scheme (SaSID, SaS2D). 



6.2 Parton parametrizations for the virtual photon 

The Q 2 -dependence of the parton distributions of the virtual photon (P 2 7^ 0) for 
A 2 <C P 2 -C Q 2 follows from the corresponding evolution equations 0], as for the 
real photon (P 2 = 0) case. Experimentally important is, however, the low-P 2 region 
A 2 < P 2 . The parametrizations g|f, ||, valid for < P 2 , are constructed in such 
a way that they reproduce the distributions of the real photon in the limit P 2 — > 
and obey the exact Q 2 -evolution equation in the region A 2 C P 2 . 
Schuler - Sjostrand 

An extention of the SaS parton distributions in the real photon to the virtual photon 
case. In the point-like contribution, the integral over the virtuality of 7* — > qq state, 
k 2 , is modified by a factor ( fc2 ^ 2 p2 ) 2 . In the hadronic contribution a factor { ^prp^ ) 2 is 
introduced, where my is a vector-meson mass. 
Drees - Godbole [fjjf 

Sets of the parton distributions obtained from the corresponding distributions for the 
real photon (DG) by including multiplicative factors to include properly the e — > ej* 
vertex. 

Gliick - Reya - Stratmann (GRS) 

The LO and NLO distributions obtained by solving the Q 2 -evolution equation with 
the boundary conditions being a smooth interpolation betwen the boundary condi- 
tions valid at P 2 = and for P 2 ^> A 2 . The applicalibility is assumed for the ranges: 
P 2 < 10 GeV 2 , 10" 4 < x < 1, 0.6< Q 2 < 5 ■ 10 4 GeV 2 and Q 2 > 5P 2 . 
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